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Summary 



In recent years the role of space weather forecasting has grown tremen- 
dously as our society increasingly relies on satellite dependent tech- 
nologies. As solar activity is the foremost important aspect of space 
weather, the forecasting of flare and CME related transient geomag- 
netic storms has become a primary initiative. Minor magnetic storms 
caused by coronal holes (CHs) have also proven to be of high impor- 
tance due to their long lasting and recurrent geomagnetic effects. In 
order to forecast CH related geomagnetic storms, the author devel- 
oped an automated CH detection package (Coronal Hole Automated 
Recognition and Monitoring; CHARM) to replace the user-dependent 
CH boundary intensity thresholding methods used in previous stud- 
ies. CHARM uses a local intensity thresholding method to identify 
low intensity regions in extreme ultraviolet or X-ray images. CHs are 
known to be regions of "open" magnetic field and predominant po- 
larity, which allowed the differentiation of CHs from other low inten- 
sity regions using magnetograms. An additional algorithm (Coronal 
Hole Evolution; CHE VOL) was developed and used in conjunction 
with CHARM to study individual CHs by tracking their boundary 
evolution. It is widely accepted that the short term changes in CH 
boundaries are due to the interchange reconnection between the CH 
open field hues and small loops. In order to test the interchange 
reconnection model, the magnetic reconnection rate and the diffu- 
sion coefficient at CH boundaries were determined using observed CH 
boundary displacement velocities. The results were found to be in 
agreement with those determined by the theory. Our results also 
indicate that the short-term CH boundary evolution is caused by ran- 
dom granular motions bringing open and small closed field lines close. 



and thus accommodating the magnetic reconnection and displacement 
of field lines. The Minor Storm (MIST) algorithm was developed by 
the author to build on the CHARM package, providing a fast and 
consistent way to link CHs to high-speed solar wind (HSSW) periods 
detected at Earth. This allowed us to carry out a long-term analysis 
(2000-2009) to study the relationship between CHs, the correspond- 
ing HSSW properties, and the geomagnetic indices. The study found 
a strong correlation between the velocity and proton plasma tem- 
perature of the HSSW stream. This indicates that the heating and 
acceleration of the solar wind plasma in CHs is closely related, and 
perhaps are caused by the same mechanism. Many authors accept 
this mechanism to be Alfven wave heating. The research presented 
in this thesis includes the small scale analysis of individual CHs on 
time scales of days, which is complemented with large scale analysis 
of CH groups on time scales of years. This allowed us to further our 
understanding of CH evolution as a whole. 
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Chapter 1 



Introduction 



This Chapter introduces the Sun - its interior structure, and the solar atmosphere 
that eoctends into the vast expanse of interplanetary space. The outer solar atmo- 
sphere, the corona, is described using a hydrostatic model in conjunction with the 
Heliospheric model of the solar wind. The origins of the solar wind are discussed 
with a more in-depth overview of coronal holes, the sources of high-speed solar 
wind and geomagnetic disturbances. Since coronal holes are the fundamental sub- 
ject of the present thesis, their observation, morphology and physical properties 
are presented together with models describing their evolution. The Chapter ends 
with the overview of space weather and the sources of geomagnetic storms. 
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1.1 The structure of the Sun 
1.1.1 The solar interior 

Our Sun is a main sequence star with a spectral type G2V, classified as a dwarf 
star with 5200-6000 K surface temperature and a yellow visible colour. Like 
other stars, the Sun is made out of high temperature mostly ionised plasma. At 
its present age of ~4.5xl0^ years the Sun still holds most of its fuel reserves 
consisting of 90% hydrogen and ~10% helium, with traces of other heavier ele- 



ments (Priest , 1984). It might make one wonder how this extraordinary hydrogen 
burning furnace we call our Sun manages to stay stable instead of going out in an 
almighty explosion. A very simplified explanation lies in the equilibrium between 
the solar gas pressure and gravity, which is also responsible for the stratification 
of the solar interior. This stratification and hydrostatic equilibrium is what en- 
sures that the solar fuel is burned gradually over thousands of millions of years. 
The temperature in the solar core is ~ 1.5x10^ K which is hot enough to ignite 
thermonuclear reactions which produce 99% of the total solar energy output. 
Nevertheless, it takes milions of years for this energy to be transferred out to the 
solar surface. In different layers of the solar interior different means of energy 
transfer processes dominate. The layer above the core is the highly opaque ra- 
diative zone (0.25-0.68 R©) where the energy is constantly absorbed and emitted 
by photons, thus the transfer happens mainly by radiative diffusion. The rigid- 
rotating radiative zone is separated from the above lying differentially rotating 
layer by a largely sheared region called the tachocline. 

From 0.68-1 Rq lies the convection zone where the energy is transferred by 



2 



1.1 The structure of the Sun 



Thi rd tier (granules) ■ 
Second tj«r (supergranules) 
CONVECTION ZO^ 
First tier . 



' Ultraviolet rays - 



RADIATION ZME : 




Random iMal1( 
radiatton pattern i 



.PHOTOSPHERE 



Figure 1.1: Schematic of the solar interior with structures of the solar atmosphere. 



kinetic motions. Here, the temperature gradient is so high that the plasma cannot 
remain in static equihbrium and convective instabihty occurs aUowing the plasma 
to rise and fall in a cellular circulation, creating the structures better known as 



supergranular cells (see Figure 1.1) 
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1.1.2 The solar atmosphere 
1.1.2.1 The photosphere 

At approximately 1 R© is the commonly accepted boundary between the solar 
interior and atmosphere. This region is the solar photosphere, also popularly re- 
ferred to as the 'surface' of the Sun since historically, it was the first part of the 
solar atmosphere to be observed with the most basic techniques (i.e., projection, 
or direct observation using a brightness decreasing filter). The observations show- 
ing the photosphere are the so-called white-light images that reveal sunspots (also 
known as active regions in corona images), which are the sites of solar magnetic 
activity. The magnetic field in the quiet Sun is ±5 — 10 G while in sunspots it can 
be well over 1000 G. The average temperature of the photosphere is ~ 5800 K 
while sunspots are relatively cooler with temperatures ~ 4500 K. This lower tem- 
perature also means a lower intensity in sunspots and makes them appear darker 
than the rest of the Sun. Unlike sunspots, faculae are brighter than the rest of the 



solar photosphere (Figure 1.2). This can be explained by reduced plasma density 
due to the strong magnetic fields conglomerated at the granular lanes. The low 
density plasma is nearly transparent, and hence allows us to see deeper layers, 
where the gas is hotter and brighter. At high spatial resolutions the granulation 
pattern caused by plasma convection can also be observed. In these features the 
brighter, middle part of a granule is the rising hot material, while the darker sides 
are the downward moving cooler plasma. 

The majority of solar radiation comes from the photosphere creating a con- 
tinuous spectrum with dark absorption lines, due to the absorption of radiation 
in the atmospheric layers above. Most of the absorption lines are formed in the 
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Figure 1.2: Sunspots and faculae observed by the Swedish solar Telescope. 

upper photosphere, but one of the most prominent hues - the hydrogen a (Ha) - is 
created in the lower chromosphere, and can be used to determine the temperature, 
magnetic field strength (based on the Zeemann-splitting of the absorption line) 
and plasma motions in the line-of-sight direction (using the Doppler broadening 
of the line). 

1.1.2.2 The chromosphere 

The chromosphere was named after its colorful (mainly pink, due to the dominant 
red Ha line) appearance during total solar eclipses. From the bottom to the top 
of the chromosphere the temperature rises from ~ 6000 K to ~ 20 000 K and the 
density drops by almost a factor of one million. This leads the magnetic field to 
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dominate over gas motions, which can also be described by a parameter known 
as the plasma beta given by the equation, 



^ Pb ^VSvr' ^ ^ 



where Pgas is the gas pressure, Pb is the magnetic pressure, rif. is the electron 
density, /c^ is the Boltzmann constant, T is the temperature and B is the magnetic 
field strength. 

Chromospheric phenomena are best observed at Ha and Ca ll K lines. The 
most prominent chromospheric features are prominences and filaments, which are 
essentially the same solar structures appearing as bright formations above the 



limb and as dark thin lines on-disk, respectively (Figure 1.3). They are magneti- 
cally bound plasma "clouds" that due to the motion and shearing of the magnetic 
field lines can build up magnetic stress to the point where magnetic reconnection 
occurs and the magnetic field restructures through a solar eruption known as a 
coronal mass ejection (CME). Further interesting chromospheric phenomena are 
the thinly distributed jets, the spicules^ which are best visible at the solar limb. 
These jets have velocities of ~25 km s^^ and they rise to an average of 11 000 km 
before the material falls back. The lifetime of spicules is in the range of minutes 



(Suematsu et al, 2008). 
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Figure 1.3: Filament eruption observed by the Solar Dynamic Observatory At- 
mospheric Imaging Assembly. 



1.1.2.3 The transition region 

The transition region is a thin and inhomogeneous atmospheric layer (~100 km) 
between the chromosphere and the corona, where the temperature rises rapidly 
from 10^ K to over a million. It was previously modelled to be static, but high- 
resolution imaging has proved it to be a highly dynamic and finely structured 



region that can be easily disturbed (Doyle et al, 2006a Gallagher et ai, 1999b). 



This is due to upward propagating waves and heated plasma fiows, as well as 
cooling downward fiows. These dynamic processes make the transition region be- 
have as an energy source and sink with regards to the atmospheric layers above 
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and below. The high temperature gradient in the transition region makes the 
plasma highly conductive which produces strong emission in the ultraviolet (UV) 



and extreme ultraviolet (EUV) range (Dowdy et al, 1986). In a simplified model 



of the transition region the plasma is confined into filamentary strands of mag- 
netic field with magnetic pressure dominating gas motions over gravity and gas 
pressure. It is important to note, however, that besides the steep temperature 
and density gradient, the modelling of the transition region is complicated due 
to the rapid time variations in the plasma /3 parameter, and the transitions from 
optically thick to optically thin regimes as well as transitions from partial to full 



ionisation (Aschwanden, 2004). 



1.1.2.4 The corona 

The corona is the layer of the solar atmosphere that can only be seen by the 
naked eye during total solar eclipses or by using coronographs, instruments that 
block out the bright solar disk allowing the much fainter corona to be visible. 
The shape of the corona changes according to the time of the solar cycle, i.e. 
the global magnetic configuration of the Sun. At solar maximum the corona 
is visible all around the limb due to the shrunken polar coronal holes and the 



wide latitudinal range at which active regions are located (see Fig. 1.4). During 
the solar minimum the active regions appear closer to the solar equator and the 
large open polar magnetic fields (polar coronal holes) cause the bright coronal 
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Figure 1.4: High-resolution image of the solar corona during a solar eclipse ob- 
served in 2006. Photo courtesy of Miloslav Druckmiiller, Brno University of Tech- 
nology. 

structures (e.g. streamers) to appear closer to the equator. The corona can also 
be observed on the solar disk at EUV or X-ray wavelengths using space-borne 
instruments. These observations reveal structures of million degree temperatures 
such as loops and active events like filament eruptions (flares, CMEs). Coronal 
holes can also be observed at these wavelengths due to their low density and 
reduced brightness compared to the other quiet regions of the Sun. 



As demonstrated in Figure 1.5, the density in the corona is at least 10^° times 



lower than in the photosphere, while its temperature is vastly larger, reaching 
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over a million degrees. If heating was only possible through thermal conduc- 
tion, the temperature of the corona would gradually drop off with distance from 
the ~5800 K solar surface (i.e. the photosphere). However, the corona reaches 
temperatures in excess of 10^ K over a short distance of ~100 km, through the 
previously mentioned transition region. In order to maintain such high tempera- 
tures the corona has to be continuously heated. Although it is not yet confirmed 
what mechanism(s) is responsible for the high temperature of the corona, there 
are several theoretical approaches. A number of models suggest that Alfven waves 
are created during chromospheric magnetic footpoint motions, which deposit their 
energy by dissipating through magnetic reconnection, current cascade, magneto- 



hydrodynamic (MHD) turbulence or phase mixing ( Cranmer et al. . 2007 Suzuki 



fc Inutsukaj [2006] [Verdini fc Vdlil [20071 |Verdini et alj [2009] [Wang fc Sh^i% 



1991). In these models the difference between slow and high-speed solar wind 
velocities are explained by differing rates of flux-tube expansion. From an ob- 
servational point of view, the large number of impulsive small-scale phenomena 
such as nanoflares and microflares have been shown to produce energy compara- 
ble to the radiative loss of the quiet Sun (QS). This is in agreement with another 
group of heating models, which suggests that the primary source of the corona 
heating lies in the small-scale magnetic reconnection processes occurring between 
open and closed field lines in the transition region and lower corona. Here, the 
difference between slow and fast solar wind is explained by the differing rates 
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Figure 1.5: The electron density (A'^e) and temperature (Tg) profile of the solar 



atmosphere (Gallagher, 1999). 



of flux emergence, reconnection and coronal heating in different regions of the 



Sun (Fisk 2003 Fisk & Zurbuchen, 2006; Schwadron & McComas 2003; Woo 



et al. . 2004). The high temperature of the corona also means that the corona 



plasma mainly consists of electrons, protons and highly ionised heavy elements 
(e.g. Fe XII, Ca XV, Mg x etc.), that continuously stream into interplanetary 
space due to the large pressure gradient as a function of distance from the Sun. 
This continuous flow of solar particles is the solar wind. 

In Sections 11.61 and 11.31 we discuss the theoretical models that describe the 
vast extent of the solar corona and the formation and heliospheric morphology of 



the solar wind. 
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1.2 Hydrostatic model of the corona 

There are many different temporal and spatial scales to solar activity (i.e. solar 
cycle, ffares and CMEs), nevertheless the solar corona can be approximated as 
stable on a large scale which allows us to understand how the corona can extend 
into the heliosphere as predicted by Chapman and Parker in the 1950s. 

In 1957 the basic model of a static solar corona was described by [Chapman 



& Zirin (1957) by assuming a static atmosphere where the energy transfer was 
purely by conduction. For simplicity, the lower boundary of the corona was set 
to be at ro — IRq above the solar surface. As the corona was suggested to be 
highly thermally conductive, and hence the equation 



5 



K = kqT,-^ (1.2) 



was used to describe thermal conduction (Spitzer, 1962), where Tg = 10 K is 
the plasma electron temperature and = 8 x 10^'' W K^^ m~^ is the thermal 
conductivity at the lower boundary. As there is no particle flow in the assumed 
hydrostatic equilibrium, the temperature gradient yields heat flux given by 



dT 

g = -.-, (1.3) 
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where we also require Vg = 0, as we assume that there are no energy sinks or 
sources. By assuming spherical symmetry, the divergence of the heat flux yields 



We require the coronal temperature to tend to zero at large distances, so the 
boundary condition will take the form 



Te{r) -ToCff, (1.5) 



where ro ~ 7 x 10^ km is the lower boundary at the base of the corona. Chapman 
assumed the corona to be in hydrostatic equilibrium described by 



dp GMq 



where p is the pressure, G — 6.67 x 10~^^ m^ s~^ kg~^ is the gravitational constant, 
Mq — 2 X 10^° kg is the mass of the Sun, and p is the plasma density that can 
be given as 

p — n{me + rUp) ~ nrUp, (1.7) 

where p is the mass density, n is the particle number density and rUe and rUp 
are the electron and proton mass. For simplicity a pure hydrogen atmosphere is 
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assumed and thus rrie is neglected. Assuming that the electrons and protons in 
the plasma have the same temperature, the coronal pressure can be written as 



P 



2nkBT 



(1- 



where n ~ 2 X 10^^ m ^, A;5=1.38x 10 ^'^ m^ kg s ^ is the Boltzmann constant 
and /i is the mean atomic weight (equaling unity for pure hydrogen). As a final 



step Te, p andp are substituted into the pressure-balance equation (Equation 1.6) 
which gives 



p(r) = po exp 



7 GMqTUp 
5 2/cBToro 



5 

ro\7 



[1.9) 



Note that as r goes to infinity, the coronal pressure tends to a finite constant 
value {poo = 10-6 - 10-1° Pa; |Golub fc Pasachoff|1997D . Although this solution 
demonstrates the vast extent of the corona, it fails by assuming the solar corona 
is at hydrostatic equilibrium at large distances. The interstellar gas pressure 



cannot contain the finite pressure given by the above equation. It was Parker 



(1958), who pointed out a year later that stellar corona plasmas must go through 



continuous hydrodynamic expansion to allow for the pressure to drop off with 



distance and balance with the interstellar gas pressure (see Section 1.3.1). 



14 



1.3 Solar wind 



1.3 Solar wind 

The solar wind is a continuous stream of electrons, protons and ionised particles 
originating from the Sun. solar wind flows emanating from QS regions travel at 
speeds of 300-400 km s~^, while streams originating in coronal holes reach speeds 
up to 800 km s~^. Transient solar eruptions such as flares and CMEs can lead to 
plasma velocities over 2000 km s~^. 

The existence of the continuous particle stream from the Sun was first sug- 



gested by Biermann (1951) based on his study on the direction of comet tails. 



After Chapman's model of the extended corona, it was Parker (1958) who gave 



the first physical model to describe the particle outflow from the Sun. 
1.3.1 Parker solar wind model 

Due to the high density in the solar interior all plasma motions are dominated by 
gas pressure. In the solar corona the density is at least 10^° times lower and thus 
magnetic fields dominate all gas motions. However, at distances over 10^ Mm the 
plasma P once again becomes larger than one, and hence plasma pressure becomes 



dominant over magnetic pressure (see Figure 1.6). Since the magnetic field lines 



are "frozen-in" to the plasma (as explained later, in Section 1.3.2), the expanding 



plasma carries the magnetic field into interplanetary space and the rotation of 



15 



1. INTRODUCTION 



^ 10 - 



V 



— I — ■ ■ 1 1 >IIPI 




10 ^ 10 ^ 10"' 

Betati^nnKT/e^) 

Figure 1.6: The variation of the plasma /? parameter as a function of distance in 



the solar atmosphere (Gary, 2001). 



the Sun subsequently deforms the magnetic field lines into an Archimedean spiral 



or Parker spiral (Section 1.3.2) 



In the Parker solar wind model the atmosphere continuously expands and 
therefore we start with the equation of mass conservation (V ■ (pu) = 0), written 
assuming a spherically symmetric radial outfiow: 



1 d 

y.2 f^j. 



{r'^pu) = 0. 



;i.io) 



Next, we consider the momentum equation {p{u ■ V)n = — Vp + pg, with g 
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—GM(7)/r^), which can be written as 



du dp GMq 
dr dr ^ 



p^-.=-±-p—r^^ (1-11) 



where u is the radial solar wind velocity and the first and second term on the right- 
hand side refer to momentum exerted by gas pressure and gravity, respectively. 
To relate the pressure and density, we use the ideal gas law: 

p = pRT. (1.12) 

Here, we consider an isothermal corona by setting T to be constant. Using the 
previous equations and the isothermal sound speed {uc = (p/pY^'^), we deduce 
the equation of motion. 



Uc'^\ du , ^ , , 



where the critical distance (r^.) is defined as 



r=^ (114) 



In Eq. 1.13 if r = Tc then either du/dr = or -u = tic, and \i u = Uc then 



either du/dr = oo or r = Tc. From this it is clear that the critical distance 
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is of great importance, since it is where the plasma velocity reaches the sound 
speed. At distances smaller than the critical distance the solar wind flow velocity 
remains subsonic, while at distances larger than the critical distance the flow is 
supersonic. Substituting typical lower corona temperatures (T = 1 — 2 x 10^ K) 
into Uc = (p/pY^'^ = (RT)^/^ leads to a critical velocity of 90-130 km s^^. This 
leads to a critical distance of ~5.7-11.5 R©. Note the dependence of the critical 
values on the corona temperature. 



To get the Parker equation we integrate Equation 1.13 



" ^log(^Y = 4log(-]+4'^ + C (1.15) 
Ur/ \UcJ \rrJ r 



where C is a constant. There are flve classes of solutions to the Parker equation 



(see Figure 1.7). Class I is double valued and indicates a solar plasma that does 
not extend into the heliosphere, while solution II does not connect to the base of 
the corona. Therefore, these solutions are considered unphysical. Solutions Ilia 
and Illb can be ruled out as well, since they suggest supersonic flows at the base 
of the corona, which is not acceptable as such high velocity flows would cause 
shocks in the denser lower corona. Classes IV and V are acceptable with their 
sub-sonic flows at the coronal base, but they predict very different flow velocities 
as r goes to inflnity. At large distances the Class IV solution gives a subsonic 
flow {u < Uc] aka "solar breeze solution"). Considering this solution, with r ^ oo 
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Figure 1.7: The five classes of solutions of the Parker solar wind model. 



Equation 1.15 can be reduced to the form 



log ( — ) ~ 4 log ( — , , 



;i.i6) 



which can be approximated as 



;i.i7) 



This relation and Equation 1.10 {pu = 1/r^) yields p ~ const /{tcUc), and p 
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const (Equation 1.12) for r — )■ oo. This shows that solution IV is unphysical as 



it gives a finite density and pressure at large r values, which cannot be contained 



by the extremely small interstellar pressure (~ 10 Pa; Suess, 1990). 

Assuming m 3> Mc (Class V solution), the solar wind fiow becomes supersonic 



beyond the critical point and with r — )■ oo Equation 1.15 reduces to 



4 log 



(1.18) 



which leads to 



2 [log 



1/2 



1.19) 



Equation 1.10 then yields p ~ const /(r u)^ which leads to p — ?■ for r — t- oo. 



Consequently, j9 — )■ 0, which can be matched with the low pressure interstellar 
medium. For this reason we accept the Class V solution as a physically realistic 
model. 

Finally, let us calculate the solar wind velocity at the Earth's distance from 
the Sun. We substitute r ~ 214 Rq, = 10 Rq and Uc — 90 km s~^ into the 



Parker equation (Equation 1.15) and use the Newton- Rapson method to solve 
the equation to get ~270 km s~^ as the solar wind velocity at Earth, which is in 
good agreement with the observed slow solar wind velocities. 
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1.3.2 Heliospheric structure of the solar wind 

To construct a realistic geometry of solar wind streams, we have to understand 
the relationship between the solar wind and the interplanetary magnetic field lines 
that originate from the Sun. It is well known that the interplanetary magnetic 
field lines are "frozen in" to the plasma due to high electrical conductivity. To 
demonstrate this we use the induction equation: 

dB 

— X {v X B) +7]V^B, (1.20) 

where v is the plasma velocity, r] is the magnetic diffusion, and the magnetic 
Reynolds number is defined as Rm ~ l/?7- In case — > oo (perfectly conductive 
fiuid), the magnetic diffusion is very small and therefore the induction equation 
can be written as 

— ^Vx(vxB). (1.21) 

This equation is known as the frozen-flux theorem. Under these circumstances 
magnetic diffusion becomes very slow and the evolution of the magnetic field is 
completely determined by the plasma fiow. Parker's model has shown that the 
corona continuously expands, and consequently magnetic field lines are dragged 
out by the expanding plasma. Therefore, if there was no solar rotation the solar 
magnetic field would take up a radial configuration, while solar rotation would 



21 



1. INTRODUCTION 



r sin e 




Figure 1.8: The spherical polar coordinate system with velocity components (left) 
and the plasma displacement (right). 



wind the magnetic field lines up into a spiral. Here, we consider the latter scenario. 

Let us take a spherical polar coordinate system (r, 9, (p) that rotates with the 
Sun. The velocity components of the solar wind can be written as Ur = u,ug = 0, 
and U(f) = — f^rsin^, where Q = 2.7 x 10~^ rad s~^ is the solar angular velocity, 
and 6 is defined so it is 0° at the rotation axis and 90° at the solar equator (Figure 



1.8) 



From the radial plasma dislocation geometry shown in Figure 1.8, we can 



establish that 



r sm 



dr 



Ur ' 



;i.22) 
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and substituting u^j, = —VLrsm.6 gives 

rdcf) rVL 
dr Ur 



;i.23) 



Because tlie solar wind plasma streams are equivalent to the magnetic field lines, 
the above equation describes the differential rotation of the solar wind streams in 



the heliosphere. Integrating Equation 1.23 gives the equation for the Archimedean 



spiral geometry as described by Parker 



u 

r--ro = -^(0-0o), (1.24) 



where = const for r — )■ oo, since at large distances the solar wind speed tends 



to constant value. Equation 1.24 will be later used in Chapter p\ to calculate 



the arrival of high-speed solar wind (HSSW) streams at Earth. To demonstrate 



the magnetic field configuration in the heliosphere, let us consider Figure 1.9 It 



depicts a magnetic field line deformed into a spiral arm, with a coronal hole as 



the source of the continuous solar wind stream. As shown in Figure L9i, the 
magnetic field is radial out to 2.5 Rq, and is gradually deformed with distance. 
The distance out to which all "open" field lines are taken to be radial is commonly 
referred to as the source surface, and is used as a boundary condition for magnetic 
field extrapolations to model the large scale solar magnetic field. (Note, that this 
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b). 



coronal source sphere 
, coronal hole 




Figure 1.9: Top: the interplanetary magnetic field from a poleward view, bottom: 



as viewed from the ecliptic (Lockwood et al, 1999) 



distance can be different in different models.) 

The figure shows the magnetic field line being Earth-connected, which is the 
case when the Earth passes through a high-speed solar wind stream originating 
from a CH. At Earth, the inclination of the magnetic field line relative to the 
radial direction is given by the angle ip. Note, that ip indicates the inclination 
angle of the magnetic field to the radial direction. It is 0° at the solar surface 
where the solar wind stream is radial [Bgyu = B^), and becomes larger with 



distance, as the magnetic field lines become more deformed (see Figure 1.9 1 for 



reference). Figure 1.9 d shows the same scenario from an ecliptic view, with 
the heliographic equator being tilted relative to the Earth's orbital plane (this 
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angle is known as Bq). Let us now calculate the magnetic field components 
[Bf. = B cosip, Bq = 0,Bff, = —Bsinip). We start with Gauss's law of magnetism 
(V • B = 0) written as 

i|(r^B.) = (1.25) 

Here, we get r'^Bj. = const and we can express the constant as tq^Bq so that 



Br = (1.26) 



From this we can express B^ using Equation 1.22 (where tanip = —rQsinO/ur): 



B. = -B.tanV = b/-^''-^^^. (1.27) 



Considering a single spiral arm at the equator {6 = and sin^ = 1), from the 
previous equation we get 



tan^ = -B^/Br = — , (1.28) 



where substituting Q, r=l AU= 150 x 10^ km and Ur — 300 km s ^ (the average 
slow solar wind speed), at the distance of the Earth, we get ~ 45° as the 



inclination angle of the magnetic field to the radial direction (Goossens, 2003). 



This is in good agreement with observations. (However, note that solar wind 
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Figure 1.10: The 2D structure of the Parker spirah 

originating from coronal holes are typically around 500-800 km s~^, and therefore 
the inclination angle ip would be less than the calculated 45° at Earth.) The 



spiraling slow solar wind streams are demonstrated in Figure 1.10 

The spiral geometry changes further according to the location of the current 
sheet (the surface separating the positive and negative polarity global magnetic 
field), which in an idealised case of a perfectly dipole Sun would be a projection 
of the solar equator separating the northern and southern coronal hemisphere. 
However, the position of the equatorial current sheet can be modified by active 
regions (ARs) and coronal holes (CHs) located near the solar equator forcing it 
to move to higher or lower latitudes. The interaction of low and high-speed solar 
wind streams emanating from the QS and CHs can also change the spiral struc- 
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ture. The region where the different velocity streams meet produces compressed 
plasma (the high speed plasma catches up with slow plasma and typically forms 



a shock), also known as a corotating interactive region (CIR) (see Section 1.5.1). 



The top schematic in Figure 1.11 shows the spiral geometry in 3D, with the 
magnetic field lines originating at different latitudinal heights. It also shows the 
termination shock where the solar wind is met by the interstellar medium, which 
consequently slows the solar wind down to subsonic velocities. The streaming 
interstellar medium elongates the solar magnetospheric structure (and thus the 
interplanetary magnetic field lines), similar to how the solar wind deforms the 
Earth's magnetosphere. The heliopause is the region where the solar wind pres- 
sure and the interstellar medium pressure balance out. The bottom schematic 



in Figure 1.11 shows an even larger-scale view of the interplanetary magnetic 
field (IMF). The IMF lines are shown to be further deformed and dragged out 
by the streaming interstellar medium, and the resulting "tail" is known as the 
heliosheath. However, recent results by the Cassini and IBEX satellites indicate 
that the long accepted elongated tail-hke structure of the hehosphere might be 
incorrect. It is suggested that the heliosphere is bubble-shaped with the interac- 
tion between the solar wind and the interstellar medium occurring in a narrow 



ribbon (Frisch & McComas 2010 Krimigis et ai, 2009). 



Observations by the Solar Wind Observations Over the Poles of the Sun 
(SWOOPS) instrument aboard the Ulysses satellite have shown that the ori- 
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1.3 Solar wind 



gin of the high-speed solar wind changes according to the time of the solar cycle. 



Figure 1.12 shows the solar wind velocity and the IMF polarity (blue: negative, 
red: positive) as a function of latitudinal location. The first and third images 
from the left show the high-speed velocity wind to dominate high latitudes and 
having distinctly separate polarities. 

This is in good agreement with the dipole structure of the global solar mag- 
netic field during solar minima, a time when the solar poles are dominated by 
large polar CHs (sources of high-speed solar wind). During solar maximum , the 
polarities and the slow and high-speed solar wind sources are mixed, which can 
be explained by the multipolar structure of the solar magnetic field, the large 
number of low latitude CHs, and the generally enhanced activity of the Sun. The 



diagram at the bottom of Figure 1.12 shows the sunspot number (the indicator 



of solar activity) corresponding to the SWOOPS observation periods. 
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1.4 Coronal Holes 



Due to their lower coronal temperature and density, CHs have reduced emission 
at X-ray and extreme ultraviolet wavelengths, and hence can be identified as 



dark regions (Altschuler et al. 1972; Chapman & Bromage, 2002; Del Zanna & 



Bromage, 1997 Vaiana et ai, 1976). In Figure 1.13 a CH is shown in a high 



spatial resolution image taken by the Solar Dynamic Observatory Atmospheric 
Imaging Assembly (SDO/AIA). 

Besides CH plasma properties being different from those of QS regions, their 
magnetic configuration is fundamentally different too. QS regions are dominated 
by small and large closed magnetic fields (aka magnetic carpet), while CHs are 
regions of conglomerated unipolar open magnetic field. The open field lines are 
brought up to the solar surface and are congregated by the convection occurring 
below the photosphere. When observing the origins of the CH open field lines at 
the photospheric level it is clear that they congregate at intergranular lanes as 



shown in the schematic in Figure 1.14 

As the density and pressure drops with height above the photosphere the open 
magnetic fiux tubes expand and fill out to create what is observed as a CH at 



corona heights (Cranmer, 2005 Tu et al. , 2005). Figure 1.15 shows the mentioned 



magnetic carpet as well as the open magnetic fiux concentrations expanding with 
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Figure 1.13: Coronal hole observed by the SDO/AIA instrument at 193 A. 




Figure 1.14: A schematic of the open flux structure in CHs (Cranmer, 2009). 



height to form a polar coronal hole. The CH related magnetic field concentrations 
in the photosphere have a predominant polarity, a characteristic often used to 
classify low-intensity regions as CHs. 



It is debated whether the terminology of "open" field is correct, since the field 
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Figure 1.15: The closed magnetic field of the QS magnetic carpet and the ex- 



panding flux tubes of CH regions (Tu et al, 2005). 



lines do connect back to the Sun despite their far reach into the heliosphere. Also, 



CHs are not entirely made up of open field lines (see Fig. 1.15), at the transition 



region level lots of small loops (bipoles) can be observed (Wiegelmann et al 



2005). These bipoles might have emerged inside the CH or have been admitted 



through magnetic reconnection occurring at the boundary of CHs. Whether the 
bipoles remain within the CH or leave the CH over time due to solar rotation is 
still unclear and requires a more detailed analysis (for more detail see Chapter |4]). 



Note the bright structures (bipoles) visible within the CH in Figure 1.13 



As previously discussed in Section 1.3.1 the corona plasma continuously ex 
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pands into the heliosphere and the magnetic field fines foUow tfie plasma, as 
explained by tfie frozen flux tfieorem. Witfi CHs containing predominantly open 
field lines tfiere is no constraint in tfieir extension into tfie fieliospfiere. Due to 
tfie different solar wind acceleration mecfianism (e.g. Alfven wave acceleration 



discussed later in Section 1.4.4), tfie CHs solar wind velocity is fiigfier tfian tfiat 



of tfie QS solar wind. Tfie fiigfi-speed solar wind velocities range between 500-800 



km s"^ (Fujiki et al. ^ 2005 Obridko et al 2009 Scfiwadron et al. , 2005; Vrsnak 



et al. , 2007a). Tfie composition of tfie CH related HSSW is also very different from 



tfiat of tfie slow solar wind and can be easily identified (see Section 1.4.3). Wfiile 
botfi tfie slow and fiigfi-speed streams take up tfie spiral structure in tfie fielio- 
spfiere, tfie fiigfier velocity stream is more radial and tfius tfie fiigfi-speed plasma 



runs into tfie slower stream creating a sfiock and a CIR (see Section 1.25). Upon 



reacfiing Eartfi, tfie CIR gives rise to magnetic disturbances (Cfioi et al.. 2009 



Fujiki et al. 2005 He et al, 2010 Scfiwadron k McComas, 2003). Since CHs are 



long-lived features witfi lifetimes from days to over a year, tfiese disturbances are 
recurrent and can be predicted (Cfiapter [s]). 

1.4.1 Observation of CHs 

Tfie first evidence for CHs was found by E.W. Maunder in 1905, wfio identified 



geomagnetic disturbances recurring every 27 days. Knowing tfie solar rotation 
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period he suggested that the disturbances could be caused by a solar region. 
In the 1930s using statistical methods on geomagnetic data Julius Bartels also 
suggested that the recurring geomagnetic disturbances were of solar origin and 
named them "M-regions". The first observations to indirectly image CHs were 



made in the 1940s, but they were only identified in the 1950s when Waldmeier 



(1957) revisited the Fe xiv 5303 A coronograph images. Waldmeier identified 
the recurrence of dark regions near the solar limb and derived the shape of the 
region from limb-width measurements. Based on their low intensity and round 
shape, he named the regions "locher" (the German word for "holes"). The first 
on-disk observations of CHs were made in the 1960s and 1970s by X-ray and EUV 
sounding rockets and orbiting telescopes, and were confirmed to be low-emission 



regions. Using OSO-4 observations, Munro & Withbroe (1972) found CHs to 



also have a reduced temperature and density compared to the QS. The launch 
of Skylab in 1973 brought higher quality and quantity soft X-ray observations 



(Huber et ai, 1974), and CHs were discovered to be exempt from the differential 
rotation of the underlying photosphere. CHs can also be observed in He I 10830 A 
absorption line images as bright, blurry patches. In the QS, due to the radiation 
from the overlying corona, the He I 10830 A transition levels become overpop- 
ulated, which leads to increased absorption of radiation and hence QS regions 
appear darker. In contrast, since CHs have reduced density and hence coronal 
emission, the He I 10830 A transition levels are less populated, which leads to less 
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Figure 1.16: A coronal hole visible in four different wavelengths on 8 December 
2000: a) 20 A (Yohkoh), b) 195 A (SOHO/EIT), c) 284 A (SOHO/EIT), d) He i 
10830 A (KPVT). 



absorption, making CHs appear brighter than the QS (Andretta & Jones, 1997 



Zirin, 1975). Figure 1.16 shows a CH observed on 8 December 2000 at four differ- 



ent wavelengths: at 20 A using the Yohkoh Soft X-ray Telescope (Yohkoh/SXT), 
at 195 A and 284 A using the SOlar and Heliospheric Observatory Extreme ul- 
traviolet Imaging Telescope (SOHO/EIT), and at 10830 A using the Kitt Peak 
Vacuum Telescope (KPVT). The quality of corona observations improved sig- 



36 



1.4 Coronal Holes 



nificantly with the launch of SOHO in 1995, providing 12 minute cadence and 
2.6 arcsec/pixel resolution EUV images and magnetograms. A further advance 
was made in observations with the Hinode X-ray telescope (Hinode/XRT), pro- 
viding 1 arcsec/pixel resolution and 200 s cadence. EUV observations improved 
further with the Solar TErrestrial RElations Observatory Extreme Ultraviolet 
Imager ( STEREO/ EU VI), which has a 1.6 arcsec pixel resolution and 3 minute 
cadence. The latter two satellites were both launched in 2006. In 2010 SDO was 
launched providing the highest quality EUV images to date with a 1 arcsec/pixel 
spatial, and 10 second temporal resolution. All of the above mentioned instru- 
ments allow the visual detection of low intensity regions. However, these regions 
can only be categorised as CHs once a dominant magnetic polarity is found. 

In the present thesis data from the SOHO/EIT and the Michelson Doppler 
Imager (MDI) were used as the former instrument allows the detection of low 
intensity regions, while the latter can be used to detect a predominant polarity 
within a region of interest. Furthermore, these instruments have been operational 
since 1995, providing a remarkable database spanning over a whole solar cycle, 
which is fundamental in the long-term CH study presented in the thesis. 
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1.4.2 Coronal hole morphology 

The detection of CHs in soft X-ray images made by the Skylab Apollo Telescope 
Mount (ATM) and the Yohkoh Soft X-ray Telescope (SXT) was straightforward 
since they appeared as significantly darker regions in comparison to the QS. How- 
ever, the low spatial resolution and the broad temperature range often lead to 
the boundary appearing "blended" into the QS, making its identification diffi- 
cult. The EUV observations by SOHO/EIT, STEREO/EUVI, SDO/AIA and 
the X-ray observations by Hinode/XRT have higher spatial resolution and thus 
provide a better CH boundary contrast. It is to be noted though, that due to the 
broad temperature range of the EUV and X-ray observations a large number of 
high temperature loops and other hot material become visible and can still cause 
line-of-sight obscuration effects at CH boundaries. This effect is reduced when 
observing CHs close to the solar disk centre, and increased near the limb. Also, 
during solar minimum this effect can be smaller due to the overall reduced activity 
of the Sun. He I 10830 A observations have a much narrower temperature range 
resulting in little or no obscuration effects at CH boundaries, however, the con- 
trast between QS and CH is much lower making the CH boundary identification 



more challenging (see Figure 1.16). Based on their location and lifetime, CHs 
are popularly divided into polar, equatorial and transient coronal holes. Some 
authors debate whether transient CHs is a valid group as they are essentially a 
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Figure 1.17: The "Elephant's Trunk" transequatorial coronal hole, observed on 
27 August 1996 (a), and the same region observed several rotations later on 1 Nov 
1996, with equatorial CHs (b). 

by-product of CME eruptions when plasma vacates the region (also known as 
dimming regions). However, by definition CHs are regions of conglomerated open 
magnetic fiux of similar polarity as well as reduced density, and as such, dimming 
regions qualify as transient CHs. 

The most long-lived CHs are the polar holes that start to develop near the 
time of the pole-reversal from pre-polar holes, that form at high-latitudes from 
the diffused magnetic fiux of decaying active regions. These CHs then extend to 
form large polar holes once the pole-reversal is complete. At the start of their 
growing phase polar holes tend to be asymmetric and later develop equatorial 
extensions (e.g. the well known "elephant's trunk", see Figure [lIT^ ). Polar 



coronal holes areas reach maximum during the solar minimum ( Harvey & Recely 
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2002). Equatorial CHs (Figure 1.17 d) are present throughout the solar cycle but 



are more frequent and extensive during the solar minimum when the number of 



ARs is lower (Wang et ai, 1996). 



1.4.3 Plasma properties 



Apart from low EUV and X-ray emission and a dominant polarity, CH regions can 



also be detected based on their elemental abundances (Feldman, 1998). Unlike 



QS and ARs, CHs exhibit similar abundances in the photosphere, upper chromo- 
sphere, transition region and the lower corona, while QS regions and ARs show 
significant enhancements in the abundance of low first ionisation potential (FIP) 
elements. It is commonly accepted that in the upper chromosphere low-FIP el- 
ements become ionised and high-FIP elements remain more neutral. Since CHs 
are regions with temperatures typically lower than the QS and are dominated 
by open magnetic field, the FIP abundances found near the photosphere remain 
unchanged at higher atmospheric layers. The differing abundances of QS related 
slow solar wind and CH related HSSW can be detected at 1 AU when the HSSW 



streams sweep past Earth (von Steiger et ai, 1995 Zurbuchen et ai, 2002). It 



is not yet agreed upon what processes cause this preferential ionisation in the 
chromosphere. 



Figure 1.18 shows a semi-empirical model of the electron density (solid black 
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Figure 1.18: The empirically acquired densities in the chromosphere, transition 
region, and low corona. The solid black line shows the electron number density, 
while the dot-dashed black line corresponds to the total hydrogen number den- 



sity (Avrett & Loeser, 2008). Visible- light electron number densities for CHs and 



streamers (blue and red curves, respectively) are shown for comparison (Cranmer 



2009) 



line) and total hydrogen number density (dot-dashed black curve) in a polar coro- 
nal hole (PCH) and streamers at the atmospheric heights of the chromosphere, 



transition region, and low corona (Avrett & Loeser, 2008). The blue curves are 



visible-light electron CH density measurements (Cranmer et al. , 1999; Doyle et al. 



1999 Fisher & Guhathakurta 1995 Guhathakurta et al.. 1999a). The red curves 



correspond to visible-light measurements of equatorial helmet streamers ( Gibson 



et al.. 1999 Sittler & Guhathakurta, 1999). It must be noted that although the 



streamers were measured to be about a factor of ten denser, a factor of 2-3 vari- 



ation is common due to the different number of polar plumes projected in the 
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Figure 1.19: The north polar coronal hole plasma density as a function of the 



radial distance ( Guhathakurta Sz Fisher, 1998). The black and blue lines represent 



density estimates for the northern PCH in 1995 and 1993, respectively. The crosses 



resent observations by Allen ( 1973 ) 



correspond to measurements by (Munro & Jackson, 1977), while the asterisks rep 



line-of-sight measurements (Cranmer, 2009). In Figure 1.19 the black and blue 



lines correspond to density estimates of the northern PCH in 1995 and 1993, re- 



spectively. The crosses correspond to measurements by (Munro & Jackson, 1977) 



and the asterisks represent observations by (Allen, 1973). 



The empirical and model temperatures in PCHs and fast wind streams are 



shown in Figure 1.20 The mean plasma temperatures from semi-empirical model 



of Avrett & Loeser (2008) are demonstrated with the dashed black curve, while 



the turbulence-driven coronal heating model of Cranmer et al. (2007) are shown 
with the solid black curve. The dark and light blue bars correspond to electron 
temperatures measured off-limb by the SOHO Solar Ultraviolet Measurements of 



Emitted Radiation (SUMER) instrument by Wilhelm (2006) and Landi (2008) 
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Figure 1.20: The CH plasma temperatures determined using a semi-empirical 
model (dashed line) and a turbulence-driven coronal heating model (solid line). 
The blue symbols correspond to electron temperatures (SOHO/SUMER), the red 
symbols correspond to proton temperatures (SOHO/UVCS), and the green symbols 
correspond to O'*'^ ion temperature measurements (Cranmer 2009). 



respectively. The proton temperatures were assembled by Cranmer (2004) using 
the SOHO Ultraviolet Coronagraph Spectrometer (UVCS) instrument, and the 
O"*"^ ion temperatures (based on the O vi 1032 and 1037 A spectral lines) were 
determined by Landi & Cranmer (2009) and Cranmer et al. ( 2008[ ) (open and 



filled green circles, respectively). It is clear from Figure 1.20 that the electron 
and proton temperatures are significantly lower than the heavy ion temperatures. 

Due to the high temperature of the corona, the hydrogen atom loses its elec- 
tron and exist most of the time as a free proton. Hence, the measured properties 
of neutral hydrogen are considered as valid approximations of proton properties 



below 3 R© (Allen et al., 2000). Observations by the UVCS instrument indicate 
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Figure 1.21: Slow and fast solar wind stream speeds are shown as a function of 
radial distance. Polar coronal hole proton and O"*"^ ions velocities are shown in 
green. Theoretical models of the polar and equatorial solar wind are demonstrated 



with the black line ( Cranmer et al. , 2007 ) and the open circles correspond to the ve- 



locity measurements of blobs above equatorial streamers using the SOHO/LASCO 



instrument (Cranmer, 2009). 



that in comparison to protons, O^^ ions are heated much more strongly, reaching 
temperatures similar or higher than the core of the Sun ( 2x10^ K). These ex- 
treme temperatures of the heavy ions measured in CHs can be explained by the 
preferential heating of different elements by processes, such as reconnection and 
turbulence driven heating, Alfven wave heating, and heavy ion velocity filtration 



(for details see Cranmer 2009). 



Figure 1.21 shows the solar wind speeds measured in a PCH using the UVCS 



instrument, together with the velocities derived from a theoretical model of the 



fast solar wind (Cranmer et a/., 2007). Polar coronal hole proton velocities are 



shown in red (Kohl et a/., 2006) and O^^ ions velocities are shown in green 
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Temperatures and Densities in Coronal Holes 



Data Type 


Radial Span (RJ 
of Observations 


T. 
(10* K) 


(10« cm"') 


Reference 


Limb 


White liglit 


1.04-1.4 
1.1-1.6 


1.15 
1.2 


2 (at 1.2 RJ 
2.3 (at 1.1 RJ 


Van de Hulst 1950a, b 
Saito 1965 


X-rays 


1-1.14 
1.04-1.1 


1.3 

0.77-0.81 


27 

7-8 (at 1.04 RJ 


Krieger et al. 1973 
Ahmad & Webb 1978 


XUV 


1.025-1.35 
1.1-1.3 


0.7-1.2 

0.9-1.37 (at 1.15 R, in SCH) 


2 (at 1.1 RJ 

0.62 + 0.11 (at 1.15 RJ 


Bohlin et al. 1975 
Guhathakurta et al. 1992 


EUV 


1.0-1.1 
1.0-1.15 


1.0-1.28 (at 1.15 R, in NCH) 
1.1 (at 1.06 R,) 
0.76 (at 1.03 R,) 
1.1 (at 1.1 RJ 


0.45 ±0.11 (at 1.15 RJ 
5 (at 1.03 RJ 
1.7 (at 1.03 RJ 
0.41 (at 1.1 RJ 


Ahmad & Withbroe 1977 
Mariska 1978 


UV 


1.01-1.02 


<1.0 




Doschek & Feldman 1977 


Lya 

(res. scat.) 


1.5-3.5 
1.5-2.0 
1.5-3.5 


1.22 + 0.1 (at 1.5 RJ 
1.7 + 0.3 (at 1.8 RJ 
1.34 (at 1.5 RJ 


0,07 (at 1.5 RJ 
0.5 (at 1.5 RJ 
0.017 (at 1.5 RJ 


Withbroe et al. 1985 
Withbroe et al. 1986 
Strachan 1990 


Disic 


EUV 




1.05 (at center of CH) 
1.23-1.44 (at borders of CH) 

1.1- 1.3 (in network) 

1.2- 1.6 (in cell) 


3 (at center of CH) 


Munro & Withbroe 1972 
Raymond & Doyle 1981a 


In Situ 


Charge state 


<1.5> 


1.3 ± 0.2 (CNO ions) 




Galvin et al. 1984 



<1.5> 
<3.0> 



1.3 ± 0.3 (CNO ions) 

1.4 ± 0.2 (Fe ions) 



Ipavich et al. 1986 



Figure 1.22: List of CH density and temperature measurements using on- and 



off-disk observations (Habbal et al, 1993). 



(Cranmer et al, 2008). For comparison, the observed and theoretically derived 
velocities for slow solar wind streams (associated with equatorial helmet stream- 
ers) are shown at the solar minimum. For the latter, "blobs" were tracked above 
helmet streamers using the SOHO Large Angle and Spectrometric Coronagraph 



(LASCO) instrument (Sheeley et al. 1997) 



For an extensive list of CH density and temperature measurements see Habbal 



et al. (1993) and references therein (Figure 1.22). 
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1.4.4 Solar wind acceleration theory 



In order to recreate the measured properties of CHs and the corresponding 
HSSWs, a MHD model needs to be constructed, including an appropriate coro- 



nal heating source (Tziotziou et ai, 1998). We start with the hydrodynamic 



model that includes the continuity equations of mass (Equation 1.10), momen 



tum (Equation 1.11), and energy. The latter can be written as: 



1 d 



yv{e + p)) = -pvg + Qtot + vF^ 



Alfveni 



;i.29) 



where FAifven is the force due to Alfven waves, Qtot is the sum of all external 
energy sources: 

Qtot Qop ~\~ Qrad ~\~ Qcond ~\~ Qmechi (1.30) 

with Qop being the opacity heating of the gas due to the absorption of pho- 
tospheric radiation, Qrad representing the radiation losses, Qcond is the heating 
through thermal conduction and Qmech is the mechanical heating. The e in Equa- 



tion 1.29 is the total energy density defined as 



^ 2 , P 

e= -pv + -, 

2 7 — 1 



;i.3i) 



where 7 is the ratio of specific heats (equals 5/3). We also include the ideal gas 



law (Eq. |1.12 ). Additionally, it is assumed that Alfven waves propagate along the 
magnetic field lines and exert a force on the gas flow. This force is the gradient 
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of the Alfven wave pressure. It is assumed that the Alfven waves propagate in a 
magnetic field that varies with distance as 



B = Bn 



;i.32) 



where Bq is the magnetic field at the surface of the Sun. The Alfven wave force 
can then be written as 



Alfven 



-VP, 



1 de 



;i.33) 



where is defined as the wave stress tensor and e as the energy density of the 



Alfven waves. For further details see Tziotziou et al. ( 1998 ) and references therein. 



This model successfully reproduces the high speed velocities (~630 km s^^) typ- 
ical of CH solar wind streams observed at Earth, while a simple hydrodynamic 
model would yield a velocity of ~5 km s~^. For a comparison between the proper- 
ties deduced from a model excluding and including the Alfven wave acceleration, 



together with the observed properties, see Figure 1.23 



Regarding the interaction of open and closed magnetic field lines at the bound- 



aries and inside CHs, one of the main theories is the interchange model (Fisk 



2005 Lionello et al. . 2006 Wang et al.. 1996). It is based on the assumption that 



the dominant process in coronal open field evolution is the reconnection between 
open and closed magnetic field lines (described in detail in Chapter |4]). This view 
was initiated by observations showing streamer evolution (infiows and outfiows) 
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Quantity 



Base mechanical flux Fo 
Mechanical flux Fti [a] 
Alfven base flux Fao 
Magnetic field Bo at 1 Rq 
Pressure ptr 

Maximum Temperature Tmaj 
Distance of Tma.x 
Distance of sonic point 
Wind velocity v at 1 AU 
Density p at 1 AU 
Temperature T at 1 AU 
Mass loss AI 
Emission measure EM 
Average temperature Tav [7] 



Model 
without Alfven pressure 

1.28 X 10^ ergcm^^s"^ 
2 X 10^ ergcm"^ s"^ 
none 
none 

4.6 X 10"^ dyne cm"^ 

1.9 X 10" K 

1.14 Ro 

none [5] 

5 kms^^ 

1.2 X 10"^^ grcm"^ 

2 X lO'' K 

1.9 X 10"^" Moyr"^ 
10^=-^3 ^^-5 

1.7 X 10'^ K 



Model 
with Alfven pressure 

1.28 X 10'' ergcm"^ s"^ 
2 X 10^ ergcm"^ s"^ 
1 X 10*^ ergcm"^ s"^ 
10 G 

4.2 X 10"^ dyne cm"^ 
1.7 X 10^ K 
1.14 Rq 
10 Ro 

630 kms"^ 

1.3 X 10"^^ grcm"-'* 

2 X 10'' K 

2.9 X 10"" Mgyr"^ 
lO^'^.^s 

1.5 X 10'' K 



Observations 
and previous models 



few 10'^ erg cm ^s ^ [1] 
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[a] Fti is the mechanical flux through the transition region 

Ptr is the pressure at the transition region 
[7] Tav is the averaged temperature with the emission measure as a weight function 
[5] The flow is subsonic up to 1 AU. The sonic point lies at r = 1714 Rq 
[1] Model values. See Withbroe (1988) 

[2] Inferred by Skylab observations and models. See Hundhausen (1977) 

[3] Value inferred by EUV observations and models. See Munro & Withbroe (1972) 

[4] Yohkoh observations. See Hara et al. (1994) 

[5] There are no observations. The value depends on the model considered. 
[6] Ulysses observations. See Phillips et al. (1995) 
[7] Inferred by several in situ observations. 



Figure 1.23: The CH and related HSSW properties derived from a model ex- 
cluding and including Alfven wave acceleration. The observed CH and solar wind 



properties are listed for comparison (Tziotziou et al. 1998). 



which could be explained with the closed flux opening into the solar wind and 



open flux closing back into the streamer (Sheeley & Wang, 2002). Interchange 



reconnection can also explain the rigid rotation of the CH boundaries (Timothy 



et al, 1975), by allowing closed magnetic fields rotating at the local QS differ- 



ential rotation rate to continuously reconnect with the CH open field lines and 
effectively "pass" through CHs. This would allow CHs to resist the shearing effect 
of the differential rotation at the photospheric level. At corona heights no such 
compensation mechanisms are necessary as the corona is well known to be rigidly 
rotating. The counter-argument to interchange reconnection is that the contin- 
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uous opening and closing of coronal field lines would result in the injection of 



magnetic flux into the solar wind and cause a bi-directional heat flux (Antiochos 



et ai, 2007). However, the bi-directional plasma jets observed at CH boundaries 



by Madjarska et ai (2004) does support the interchange reconnection model. 



The other main theory is the "quasi-steady model" (Antiochos et ai, 2007) in 
which the global solar magnetic field is determined by the instantaneous distri- 
bution of photospheric magnetic flux and the balance between magnetic and gas 
pressure in the corona. The disadvantage of the quasi-steady model is the assump- 
tion of a static coronal magnetic field with topologically well separated regions 
of open and closed magnetic field. Since the corona does change in time due to 
photospheric motions and flux emergence, its evolution is modelled through a se- 
ries of time stationary stages. The simplest quasi-steady model was the potential 



field source surface model (Altschuler & Newkirk, 1969), where the gas pressure 
was neglected and the magnetic field was considered current free. The magnetic 
field was then extrapolated based on photospheric flux values, assuming the field 



to be radial from the given distance of the source surface (Luhmann et ai, 2003 



Schrijver & DeRosa, 2003). In recent years, the quasi-steady model has been 



improved by including the solution of the full MHD equations, which made the 



assumption of a fixed surface and a current-free corona unnecessary (Odstrcil 



2003 Pisanko, 1997; Roussev et al.. 2003). These models even allow the relax- 



ation of the quasi-steady approach by including photospheric time-dependence. 
However, improvement is still needed since it makes the model computationally 
heavy and no flux emergence and cancellation can yet be included. 
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1.5 Space weather 

The growing dependence of our society on advanced technologies has lead to an 
increasing vulnerability to space weather effects. While our understanding of the 
physical processes connecting solar activity and geomagnetic storms has advanced 
enormously in the past fifty years, forecasting the time, duration and scale of 
Earth-bound severe space weather events still remains a challenge. First and 
foremost, the solar regions which are likely to erupt in the form of a flare and/or 
a CME have to be identified and located. Secondly, the time of eruption needs 
to be predicted, which holds probably the largest challenge. Despite numerous 
models having been developed, the success rate of flare and CME forecasts is still 
relatively low. Nevertheless, if a region is forecasted to erupt, it still has to be 
determined whether the flare or CME is Earth-bound and whether it is likely 
to be geoeffective. Once a geomagnetic storm is forecasted, the estimated time 
of arrival (ETA) of the geomagnetic storm at Earth has to be provided together 
with the duration of the storm and the estimated end-time. This would allow 
public, commercial and military organisations to take precautionary measures in 
preparation of a geomagnetic storm. 

Programs that aim to develop space weather forecasting methods and provide 
space weather warnings publicly have only emerged in the past twenty years. To 
mention a few: the European Space Agency (ESA) Space Weather Web Server, 
the Solar Influences Data Analysis Center (SIDC) operated by the Royal Observa- 
tory of Belgium, the European Space Weather Portal (ESWeP), the U.S. National 
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Space Weather Program (such as the National Oceanic and Atmospheric Admin- 
istration Space Weather Prediction Center (NOAA/SWPC)), the Russian Space 
Weather Initiatives (such as ISMIRAN), the Austrahan Space Weather Agency, 
the Solar Activity Prediction Center (China), and the organisation connecting nu- 
merous space weather centres worldwide - the International Space Environment 
Service (ISES). Many of the above agencies provide space weather forecasts, real- 
time now-casts on solar and geomagnetic indices, observed halo-CMEs as well as 
mail-alerts. 

SWPC established three scales for the differentiation of space weather events 
based on their driving source. All three categories scale from 1 to 5, correspond- 



ing to "minor" and "extreme", respectively (Biesecker et al, 2008). 
1. The Radio Blackout Scale is based on X-ray flux measurements by Geosta- 
tionary Satellites (GOES). Enhanced levels of X-ray flux cause heightened ionisa- 
tion in the day-side ionosphere that can lead to sudden ionospheric disturbances 
(SIDs), which cause interruptions in telecommunications systems. Since many 
communication systems transmit signals using radio waves reflected off the iono- 
sphere, geomagnetic storms can lead to some of the radio signals being absorbed 
in the disturbed ionosphere, while some still get reflected but at an altered angle. 
This results in rapidly fluctuating signals (scintillation), phase shifts and changed 
propagation paths. High-frequency radio bands are especially sensitive to geo- 
magnetic activity, while TV and commercial radio stations are barely affected 
due to the use of transmission masts. 

2. The Solar Radiation Storm Scale is measured from the solar energetic particle 
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Figure 1.24: During a solar eruption electrons, protons and heavy ions are ejected, 
creating a hazardous environment for astronauts and electrical systems aboard 
satellites. Since electrons are much faster, they can be detected minutes before the 
ion radiation reaches Earth, providing a useful forecasting tool. Figure courtesy of 
the Southwest Research Institute. 



(SEP) flux. These energetic particles have energies 10-100 MeV, and sometimes 
in excess of 1 GeV. SEPs accelerated during solar eruptions take minutes to 
hours to reach Earth depending on their energies. The SEP events disturb the 
ionosphere which can lead to communication and navigation problems, while the 
elevated energetic proton and heavy ion flux levels can harm astronauts during 
extra- vehicular activities (EVAs) and damage electronic devices aboard satellites. 



The top cartoon in Figure 1.24 shows a solar eruption during which electrons and 
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ions are ejected from the Sun and are heading towards Earth. The first type of 
particles to arrive are the relativistic electrons, which are followed by the ions tens 
of minutes later. Hence, the electron signal can be used to predict the radiation 
hazard posed by energetic ions, allowing astronauts to seek shelter, and electri- 
cal devices on satellites to be put into safe-mode. Such ion radiation forecasts 
are made by the Comprehensive Suprathermal and Energetic Particle Analyzer 



(COSTEP) aboard SOHO (Posner et al, 2009). Note, satellites, and spacecraft 
in low-Earth orbit are protected to a certain extent from energetic particles by the 
Earth's magnetosphere. However, radiation is still higher outside the protective 
shield of the Earth's atmosphere, and during severe solar storms astronauts on 
a space station may reach or even exceed a whole year's radiation limit. During 
extreme space weather conditions even high-altitude aircraft routes have to be 
modified due to the enhanced risk of radiation exposure to flight crews. 

3. Geomagnetic Storm Scale is based on energetic solar events, during which 
large amounts of plasma and energy departs the Sun in the form of CMEs. The 
ejected plasma contains imbedded magnetic structures, also known as a magnetic 
cloud. The magnetic fields in CMEs need to be of strong southward direction 
to reconnect with the mainly northward directed magnetic fields of the day- 
side magnetosphere. This reconnection gives rise to strong induced fields and 
electric currents in the magnetosphere, ionosphere and the solid body of Earth 



(Figure 1.25). The latter can adversely affect large-scale power grids since the 
induced electric currents can enter long transmission lines and reduce the power 
by unbalancing the grid operation, and can even destroy transformers through 
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Figure 1.25: The effects of geomagnetic storms experienced in the Earth's atmo- 
sphere and on-ground. 

heating and vibration. The most well-known example to a power-grid failure due 
to a severe geomagnetic storm was in Quebec on 13 March 1989, which caused a 
power-outage for 9 hours. 

A CME can be geoeffective even if it does not contain southward magnetic 
fields. The compressed plasma between the shock and the CME can contain a 
southward directed magnetic cloud, which can be geoeffective. Another geomag- 
netic effect of a CME is the compression of the magnetosphere on the day-side 
of the Earth. This may leave the high-altitude satellites exposed to harmful 
energetic solar wind particles. 
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Figure |1.26| shows the classic models of the magnetospheric reconnection sug- 



gested by Dungey (1963). The top diagram shows the change in the Earth's 
magnetospheric topology as an incoming south-directed IMF reconnects with the 
typically north-directed magnetic field of the Earth's day-side magnetosphere. 
Reconnection first occurs at the Sun-facing side of the magnetosphere, where the 
oppositely directed magnetic field lines meet. Then the reconnected field lines 
travel with the solar wind, which extends them to become part of the magnetotail 
where reconnection occurs again where oppositely oriented magnetic field lines 
meet. During the second reconnection two separate magnetic structures are cre- 
ated - one that connects the North and South Pole of the Earth and a disconnected 
field line that is carried away by the solar wind. During the second reconnection 
solar energetic particles are channelled towards the poles of the Earth where the 
atmospheric particles get excited and create the well known phenomenon - the 



aurorae borealis (Figure 1.27). At the same time geomagnetic disturbances may 
be experienced at high latitudes based on the severity of the storm. 

The second, geomagnetically less effective scenario is shown in the bottom of 



Figure |1.26[ Here, the incoming magnetic field is north-directed, similar to the 
Earth's magnetic field. In this case no reconnection occurs at the day-side of 
the Earth's magnetosphere, however the incoming magnetic field gets deformed 
around the magnetosphere. At the tail the incoming magnetic field encounters 
the night-side closed magnetic field of the Earth which is oppositely directed 
to it and causes and double reconnection as shown in the figure. Through this 
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Figure 1.26: The reconnection process shown in case of a southward (top) and 



northward (bottom) directed IMF impacting the Earth's magnetosphere (Dungey 



1963). 



double reconnection a close field is created that becomes part of the day-side 
magnetosphere, while a disconnected field is also created which is carried away 
by the streaming solar wind. 

It is important to note that flux is added to the geomagnetic tail during the 
reconnection process involving a southward interplanetary magnetic field, while 
flux is removed from the tail in the case of an incoming northward interplanetary 



field. 
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Figure 1.27: Aurora seen in Sommaroya, northern Norway on 8 September 2010. 
Photo courtesy of Thilo Bubek. 

1.5.1 Corotating Interactive Regions 

Although on a smaller scale, CHs can also cause geomagnetic disturbances. When 
the HSSW encounters the slower solar wind, the solar wind plasma becomes 
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compressed (creating a CIR), while at the traihng part of the HSSW a rarefaction 



region is created (Figure 1.28). Most CIRs are bounded by shocks, but mostly at 
distances of 2 AU or more (in general, waves have to go through several nonlinear 
steepening phases before a shock is formed). Since the interaction region is of 
high pressure, it expands into the plasma ahead and behind at high speed creating 
two shocks on the two sides of the pressure enhancement. This phenomenon is 
called the "forward-reverse shock pair", where one shock travels towards and 



another away from the Sun (Gosling, 1996). The propagation and the pressure 



gradient of the shock waves lead to the acceleration of the slow solar wind and the 
deceleration of the high-speed solar wind. This interaction limits the steepening 
of the stream and transfers energy and momentum from the high-speed solar 
wind to the slow solar wind. At large distances the stream pressure enhancement 
gradually reduces and the compression regions merge as they expand. 



An extensive statistical study carried out by Choi et al. (2009) has shown 
that approximately half of the CH related CIRs caused geomagnetic disturbances 
with "Dst <-50 nT (-5x10^^ G). D^j index indicates geomagnetic storm severity, 
with high values e.g. 100 nT (10^'^ G) corresponding to a calm magnetosphere, 
while low values e.g. -300 nT (-3x10^'^ G) indicate a large storm. It has also 
been shown that 70% of the CHs were located between ±30° latitude, and the 
area of geoeffective CHs were found to be larger than 0.12% of the visible solar 



hemisphere. Choi et al. (2009) used KPVT CH maps and Advanced Composition 
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Figure 1.28: The CIR and rarefaction caused in the hehosphere by the interaction 



of the slow and high-speed solar wind plasma. Modified from Parker (1965). 



Explorer (ACE) and Wind solar wind data to link CHs to CIRs. 

With real-time data available, the arrival of CH related CIRs at Earth can 
also be forecasted. This is possible since CHs are recurrent and visible for long 
periods of time, thus making the corresponding geomagnetic storms easier to 
forecast than disturbances caused by flares and CMEs, where the eruption time 
is challenging to predict. Knowing the range of typical CH solar wind stream 
velocities, the path of the streaming plasma can be calculated together with its 
arrival time at Earth. In Chapter |5] we discuss how the CH detection algorithm 
presented in this thesis can be useful in predicting minor geomagnetic storms 
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based on the early detection of CHs. 

1.6 Thesis outline 

The CH detection method presented in the thesis provides a fast, subjective and 
automated method to identify CHs and improve the determination of their physi- 
cal properties throughout the solar cycle. The automated nature of the algorithm 
allows the detection of CHs on short and long timescales which can further our 
understanding of the magnetic and morphological evolution of CHs. This reflects 
on the small scale processes governing the CH boundary evolution as well as the 
mechanism of solar wind acceleration in CHs. Linking CH observations to HSSW 
streams using in-situ solar wind data provided a practical aspect to the presented 
research by enabling the real-time forecasting of minor storms and the investi- 
gation of correlation between CH and solar wind properties. The geomagnetic 
effects of CH related HSSW streams were also compared to changes observed in 
the Earth's radiation belt. Furthermore, the observations were used in conjunc- 
tion with the solar cycle to derive the role of CHs in the solar magnetic pole 
reversal mechanism. 

Chapter |2] presents the instruments that are used in this study to detect 
CHs. The most optimal instrument pair with the largest database, SOHO/EIT 
and SOHO/MDI, were fundamental for developing the CH detection algorithm 
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and the research presented in Chapter |4] and Chapter [5| Chapter [3] describes 
previous methods of CH detection as well as detaihng the detection algorithm 
developed by the author. The algorithm is further extended to study the small- 
scale evolution of CH boundaries in conjunction with the interchange reconnection 
model, as presented in Chapter |4j The physical properties of the studied CHs 
are determined together with the diffusion coefficient and the rate of magnetic 
reconnection occurring at CH boundaries. The large-scale and long-term changes 
in CHs are discussed in Chapter |5] where the results involving the analysis of over 
nine years of EUV data are presented. The observed properties of CHs in terms of 
the solar cycle are discussed with respect to the geomagnetic effects experienced 
at Earth over the same period. In Chapter [6] the theoretical implications of the 
short and long-term evolution of CHs is discussed, and the minor geomagnetic 
storm forecasting method is evaluated with suggestions for further development 
and improvements. 
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Chapter 2 



Instrumentation 



In this Chapter, we discuss the EUV and X-ray imagers that were used in the 
present research to detect low-intensity regions. These include the SOHO/EIT, 
STEREO/EUVI and Hinode/XRT instruments. In order to determine whether 
a low intensity region is a CH, magnetic field analysis was carried out using full- 
disk magnetograms provided by SOHO/MDI. For forecasting purposes the detected 
CHs were linked to high-speed solar wind streams using in-situ solar wind data 
from the ACE/SWEPAM and STEREO/PLASTIC solar particle monitors. 



63 



2. INSTRUMENTATION 



2.1 Introduction 

Until the 1970s, imaging the solar corona was only possible from ground-based 
observatories. Although ground based solar telescopes are still used world-wide, 
space-based telescopes are now the dominating source of data. The main problem 
ground-based observations face is the atmospheric distortions due to air turbu- 
lence. An additional distortion is caused by the optics expanding due to heat 
from the environment or the instrument itself. Although there are now advanced 
adaptive optics that correct distortions due to atmospheric turbulence within 
split seconds, airborne telescopes still have the advantage of consistently provid- 
ing high resolution images independent of weather conditions. Furthermore, EUV 
and X-ray imaging can only be done from space due to atmospheric absorption at 
these wavelengths. The disadvantage of airborne telescopes is maintenance and 
repairs, which is extremely expensive and sometimes even impossible. The size 
of the telescope is also an issue as the larger the instrument, the more it costs to 
put it into orbit. 

One of the solar features that could only be imaged from space were CHs. 
Despite this, the effects of CHs were detected long before X-ray and EUV imaging 
was carried out from space. In 1905, British astronomer E.W. Maunder found 
recurring magnetic disturbances caused by sources invisible on white-light solar 
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images ( Oliver, 1994). These were later called M-regions (M for 'magnetic') by 



Julius Bartels. The existence of CHs was first confirmed by rocket-borne X- 



ray telescopes in the 1960s (Cranmer, 2009). In 1973 Skylab was launched, which 



amongst others allowed the continuous observation of CHs. In 1991 came Yohkoh 
with the Soft X-ray Telescope and in 1995 SOHO was launched enabling the EUV 
imaging of CHs and providing data for over a decade now. In recent years Hinode, 
STEREO A & B and SDO have been launched, providing ever more impressive 
X-ray and EUV images of the Sun with the highest resolution and time cadence 
to date. 

Although magnetograms do not reveal CHs like EUV or X-ray images, they 
have proved crucial in differentiating CH from filaments and other features that 



might appear as low- intensity regions at the mentioned wavelengths (Wang et al. 



1996). CHs can also be detected indirectly by analysing the solar wind properties, 
as the particle composition of CH related solar wind streams differs from that 
originating elsewhere. 



2.2 SOHO 

SOHO is a joint mission of the European Space Agency (ESA) and National 
Aeronautics and Space Administration (NASA). The main goal of the mission is 
threefold - to study the solar interior using helioseismology techniques, to reveal 
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the mechanisms that heat the solar corona and to explain how the solar wind is 
created and what causes its acceleration. 

The spacecraft was launched from the Cape Canaveral Air Station (Florida, 
United States) on 2 December 1995. SOHO is located at the first Lagrangian 
point (LI) approximately 1.5 million kilometres from Earth (four times the dis- 
tance between the Earth and the Moon). It is always facing the Sun, providing 
uninterrupted observation. The mission had a two-year lifetime, but has now 
been in operation for almost 15 years. 



There are twelve instruments aboard SOHO (Domingo et ai, 1995): three 
spectrometers that study the plasma of the solar atmosphere (CDS, SUMER, 
UVCS), three instruments detect particles of solar, interplanetary or galactic 
origin (CELIAS, COSTEP, ERNE, SWAN), two instruments that measure the 
oscillations in the solar irradiance and gravity, and global velocity (VIRGO and 
GOLF, respectively), a coronograph that studies the off-disk solar corona (LASCO), 
an extreme-ultraviolet imaging telescope (EIT), and a magnetograph that mea- 
sures the vertical motions on the solar surface, the acoustic waves in the solar 
interior and the longitudinal magnetic field of the Sun (MDI). The instruments 



aboard SOHO are shown in Figure 2.1 



Data from the EIT and MDI instruments were used extensively in this thesis 
work and are given a detailed review in this section. 
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Figure 2.1: A schematic of SOHO and its instruments (Domingo et al., 1995) 



2.2.1 EIT 



The EIT instrument ( Delaboudiniere et a/. 1995) provides full-disk solar images 



of the transition region and the inner corona out to 1.5 solar radii above the 



limb. The telescope is a Ritchey-Chretien design (Figure 2.2). The primary 



and secondary mirrors consist of four quadrants, each of which have multilayer 
coatings that define four spectral bandpasses by interference effects. At any 
time of observation, a rotating mask allows only one quadrant of the mirrors 
to be illuminated by the Sun. The multilayers consist of alternating layers of 
molybdenum and silicon. The telescope field of view is 45 x 45 arcmin with 



2.6 arcsec pixels and a maximum spatial resolution of 5 arsecs. The detector is a 
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Figure 2.2: The schematic of the EIT instrument ( Delaboudiniere et al. 1995) 



1024 X 1024 pixel backside-illuminated charge coupled device (CCD). 

The four narrow bandpasses are centered at 171 A (Fe ix/x), 195 A (Fe xil), 
284 A (Fe xv) and 304 A (He ll) are used to study solar features visible at cor- 



responding temperatures as listed in Table 2.1 The temperature response of the 



EIT filters are shown in Figure 2.3 In the present work we use 195 A images. 



Note, that the peak temperature for the 195 A filter is at 1.6 MK, and at lower 
temperatures the signal decreases significantly. This allows the lower tempera- 
ture CHs to appear dark in comparison to the higher temperature QS and ARs. 



Figure 2.4 also shows that for temperatures of ~l-2 MK the Fe xil line intensity 



is significant at 195 A while at lower temperatures of 0.3-0.5 MK there is no con- 



tribution from Fe XII (Figure 2.5). The synthetic spectra shown in Figure 2.4 and 



2.5 were created using the CHIANTI atomic database and software (for more in- 
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Table 2.1: The EIT bandpasses, peak temperatures and observables (Delabou 



diniere et al. 



1995) 



Wavelength Ion 



Peak Temperature Observational Objective 



304 A 
171 A 

195 A 
284 A 



He II 8.0 X 10^ K 
Fe ix/x 1.3 X 10*^ K 



Fe XII 
Fe XV 



1.6 X 10*5 K 
2.0 X 10*^ K 



chromospheric network, coronal holes 
corona/transition region boundary, 
structures in coronal holes 
quiet corona outside coronal holes 
active regions 



formation see http:/ /www. ukssdc.ac.uk/solar/chianti/applications/index.html). 

In the present work we used data from the 195 A wavelength band as it pro- 
vides the best contrast between the brighter QS and the darker CHs. During the 
EIT data analysis, the data was calibrated using the standard techniques provided 



by the eit.prep algorithm available in the SolarSoft software package (Freeland 



& Handy, 1998). The calibration procedures are detailed in Section 3.3.1 



2.2.2 MDI 



The MDI instrument (Scherrer et al. , 1995 ) was primarily built to probe the solar 



interior using helioseismological techniques and measure the line-of-sight mag- 
netic field, line and continuum brightness at photospheric levels. The detector is 



a 1024 X 1024 pixel CCD, which can operate at 2 or 0.6 arcsec pixel resolution. A 
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Figure 2.3: The temperature response of the 171 A, 195 A, 284 A and 304 A EIT 



filters ( Delaboudiniere et al., 1995). 



Lyot filter and two Michelson interferometers are used to create sets of filtergrams 
obtained at five wavelengths which are combined onboard to produce the con- 
tinuum intensity maps (white-light images), velocity maps (Dopplergrams) and 
magnetograms. The line-of-sight velocity is determined from the Doppler shift 
of the Ni I 6768 A solar absorption line, while the line-of-sight magnetic field is 
determined from the circular polarisations of the v\\ component of the velocity. 
The latter allows the magnetic field intensity maps (magnetograms) to be created 
from the filtergrams based on the physical process known as the Zeeman effect. 
This is explained in the following section. 
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Figure 2.4: Spectral line intensities between 190-200 A, at IMK and 2MK (top 
and bottom, respectively). 
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Figure 2.5: Spectral line intensities between 190-200 A, at 0.3MK, 0.5MK (top 
and bottom, respectively). 
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2.2.2.1 Zeeman effect 

The effect occurs when external magnetic field is applied to a plasma, which 
results in sharp spectral lines splitting into multiple closely spaced lines due to 
the interaction between the magnetic field and the magnetic dipole moment of 



electrons. The left side of Figure |2.6| shows the single hydrogen line produced 
due to the electron transfer from the higher energy level to the lower when no 
magnetic field is present. This transition marked with a green arrow from L = 2 
to L = 1, where L is the angular momentum of the electron. The right side 
of the figure shows three lines produced due to the transitions from the energy 
levels created by the different magnetic dipole moments (m) of the electron when 
magnetic field is present. 

The number of lines produced depends on the alignment of the magnetic field 
to the line-of-sight. If the magnetic field is transverse, a line is produced at 
the original line position (vr, perpendicular to the magnetic field) and two other 
offset lines are created which are linearly polarised (cr, parallel with the magnetic 
field). If the magnetic field is longitudinal, two offset lines are produced which are 
circularly polarised (a) in opposite directions. Longitudinal magnetic fields can 
be shown by the circular polarisation of the produced a lines. The wavelength 
difference between the vr and a a spectral line is proportional to the magnetic 
field present: 

A Xb = 4.7 xlQ-^^ gX'^B. (2.1) 
Here, B is the magnetic field strength and g is the Lande factor of the spectral line. 
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Figure 2.6: A schematic of the normal Zeeman hne-sphtting. When no magnetic 
field (B) is present in the plasma, the electron transitions produce one line (shown 
on the left), while multiple lines are produced with magnetic field present (as shown 
in the right). In the normal Zeeman-effect three linearly polarised lines appear, two 
a and the middle vr when the magnetic field is perpendicular to the line-of-sight 
(LOS). Only the two a lines are seen with opposite circular polarisation when the 
magnetic field is along the LOS. 



The detection of the Zeeman splitting can be difficult or even impossible for weak 
magnetic fields. For this reason the polarisation properties are used to determine 
the magnetic field strength. The MDI instrument operates on the principle of 



the standard Babcock magnetograph (Figure 2.7). First, the a line components 



are separated using an electrooptic cell, which converts the a components into 
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Figure 2.7: The schematic of the Babcock type magnetograph operation. The 
circularly polarised light is converted into plane-polarised light by an electrooptic 
crystal, and is shifted by 90° based on the sign of the voltage applied. The Polaroid 



then selects the components (Gallagher, 1999). 



linearly polarised beams at right-angles. Next, a linear polariser is used to let 
through one beam at a time, using different sign voltages. 

The intensity difference between the signals in each state is directly propor- 
tional to the line-of-sight magnetic field strength, and hence can be used to create 
magnetograms. 



2.3 STEREO 

The STEREO A & B twin satellites were launched on the 25 October 2006 from 
Cape Canaveral Air Force Station in Florida. The mission includes two identi- 
cal satellites - the STEREO-A spacecraft that orbits ahead and slightly inside 
of Earth's orbit, and the STEREO-B that orbits behind and slightly outside the 
Earth's orbit. For this reason the STEREO satellites provide a unique opportu- 
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nity to carry out stereoscopic observations of the Sun. The spacecrafts drift away 
in opposite directions from Earth by 22 degrees per year due to the mentioned 
small differences in their distance from the Sun. When the spacecrafts are fur- 
thest from each other (being at the opposite sides of the Sun, like at the time of 
writing) they provide a nearly 360° angle view of the Sun. This has been achieved 
for the very first time in the history of solar observation. 

The STEREO mission goals are to study the initiating mechanisms of CMEs 
and their propagation through the heliosphere, to determine the location of ener- 
getic particle acceleration in the low corona, and to determine the structure of the 
solar wind. There are four instrument packages aboard STEREO, the names of 
which speak for themselves: the Sun Earth Connection Coronal and Heliospheric 
Investigation (SECCHI), the In-situ Measurements of Particles and CME Tran- 
sients (IMPACT), the Plasma and Suprathermal Ion Composition (PLASTIC) 
instrument, and the interplanetary radio disturbance tracker STEREO/WAVES 
(S WAVES). The SECCHI package contains two coronagraphs (CORl & 2), a 
heliospheric imager (HI), and an EUV imager (EUVI). 

In this thesis we used data from the EUVI and PLASTIC instruments, which 
are discussed in detail in the sections below. 



2.3.1 EUVI 



The STEREO/EUVI instrument (Wuelser et al. 2004) works on the same prin- 



ciple as the SOHO/EIT instrument in creating images at extreme ultraviolet 
wavelengths (171 A, 195 A, 284 A, and 304 A). Also a Ritchey- Chretien tele- 
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scope in its optical design, it images the Sun at the same four wavelengths as 
SOHO/EIT, but at higher image cadence and spatial resolution. 

The primary and secondary mirrors consist of four quadrants, all of which 
are coated with a narrow-band reflective coating optimised for one of each four 
wavelengths mentioned. The coating material is MoSi. Two metal film filters 
are used at both the entrance aperture and the focal plane of the telescope to 
reduce solar heat, and filter out the undesired optical, EUV and IR radiation. For 
the short wavelength (171 A, 195 A) quadrants this filter is an aluminium-on- 
polyamid foil on a coarse nickel grid, while on the longer wavelength quadrants 
a single layer of aluminium foil was used on a fine nickel mesh. The temperature 



response of the different filters is shown in Figure 2.8 The EUVI full field of view 
displays the Sun out to ±1.7 solar radii, with a 1.6 arcsec pixel resolution. The 
detector is a backside illuminated 2048 x 2048 pixel CCD. 

Although the above description refers to one instrument, two identical EUVI 
instruments were built for the STEREO mission. One is located on STEREO-A 
and the other on STEREO-B, allowing the 3D imaging of the solar corona. 

2.3.2 PLASTIC 



The Plasma and Suprathermal Ion Composition (PLASTIC; [Galvin et al. 2008) 



aboard STEREO makes in-situ measurements of the solar wind elemental compo- 
sition, ionic charge states, ion bulk flow parameters, and heliospheric suprather- 
mal ions in the energy range from ~0.3 to 80 keV/e (including ions created in 
local shock acceleration processes). These measurements were designed to inves- 
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Figure 2.8: The temperature response of the EUVI instrument (Wuelser et al. 



2004). 



tigate the link between coronal and solar wind ion processes, and the link between 
solar wind and heliospheric processes. 

PLASTIC has three fields-of-view and four instrument sections to observe 
both plasma and suprathermal particles. The PLASTIC Solar Wind Sector 
(SWS) provides a 45°x20° Sun-Earth line centered field-of-view and observes 
in the 0.3-10.6 and 0.3-80 keV energy range. The SWS consists of two parts: 
the Small Channel measures the density, velocity, thermal speed and the al- 
pha/proton ratios of the solar wind protons, while the Main Channel measures 
the composition, charge state, and the bulk and thermal speeds of the more 
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abundant solar wind minor ions (e.g., C, O, Mg, Si, and Fe). 

The two channels of the Suprathermal Ions Wide-Angle Partition Sector 
(WAP) cover a 50° x 6° and a 175° x 6° field-of-view in the off-Sun azimuth di- 
rection. Both channels measure the composition, flux and direction of suprather- 
mal ions in the 0.3-80 keV energy range. Similarly to all STEREO instruments, 
identical PLASTIC instruments were built for STEREO-A and B. 



2.4 Hinode 

The Hinode (formerly known as Solar-B) spacecraft was launched on 23 Septem- 
ber 2006 from the Uchinoura Space Center in Kagoshima, Japan. The Hinode 
mission objectives are to study the magnetic field generation and transport, the 
processes responsible for the heating and structuring of the chromosphere and 
corona, and the mechanisms triggering flares and CMEs. The three instruments 
aboard Hinode are the Solar Optical Telescope (SOT) (including a Narrow and 
Broadband Filter Instrument, and a Spectropolarimeter) , the Extreme Ultravi- 



olet Imaging Spectrometer (EIS) and the X-Ray Telescope (XRT; Golub et al. 



2007). In the present research data from the XRT instrument was used. 



2.4.1 Hinode/XRT 

The design of X-ray telescopes is radically different from classic reflectors since 
the X-ray photons are of such high-energy, they would pass through a standard 
reflector mirror. To capture the photons, they must be bounced off at very low 
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Figure 2.9: Schematic of the grazing incidence technique, which allows X-ray 
imaging. 

angles. This technique is called grazing incidence, where the mirrors are mounted 



slightly off the parallel line of the incoming X-rays (see Figure 2.9). The X-rays 
are then focused and recorded with the CCD camera. On XRT this is carried out 
with a modified Wolter I grazing incidence telescope, which has a 30 cm aperture. 

To filter out visible light and reduce heat, nine different aluminium based 
pre-filters are used made of aluminium, carbon, titanium and beryllium. The fil- 
ters provide five passbands in the 6-60 A range allowing the observation of solar 
features in the temperature range ~l-30 MK. The temperature responses of the 



different filters is shown in Figure 2.10 This shows that apart from the first two 
filters (Al-mech and Al-poly), all the filters are sensitive at very high tempera- 
tures (> 1 MK), while this sensitivity drops considerably at lower temperatures 
(~ 0.8 MK - typical of CHs), which provides a high contrast between CHs and 
the QS. The instrument has a 35 arcmin field of view with a 1 arcsec pixel spatial 
resolution. The detector is a back-illuminated 2048x2048 pixel CCD. 
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Figure 2.10: The temperature response of the XRT instrument for all the focal- 
plane filters (A = Al-mesh, B = Al-poly, C = C-poly, D = Ti-poly, E = Be-thin, 
F = Be-med, G = Al-med, H = Al-thick, and I = Be-thick). Figure courtesy of 



Golub et al. 



(2007). 



2.5 ACE 



ACE (Stone et a/., 1998) was launched on 25 August 1997 from the Kennedy 
Space Center in Florida. Similarly to SOHO, ACE orbits in the LI Lagrangian 
point observing the Sun without interruption from a distance of 148.5x10^ km 
and should be able to maintain its orbit until 2024, based on its fuel reserves. 
Due to its continuous ability to sample the solar wind, ACE can provide space 
weather warnings shortly (~1 hour) before a magnetic storm hits Earth. Often 
it is difficult to provide earlier warnings in the case of large geomagnetic storms 
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originating from flares and CMEs, liowever minor storm warnings can be improved 
by detecting the features responsible (e.g. CHs) on solar images several days in 
advance. The HSSW forecasting method discussed in this thesis relies on the data 
from the Solar Wind Electron Proton Alpha Monitor (SWEPAM) instrument 
detailed in the following section. 

ACE has three monitoring instruments and six high-resolution sensors which 
detect particles of solar and galactic origin, measure electron and ion fluxes, and 
determine the direction and the strength of the interplanetary magnetic fleld. 
These instruments are the Electron, Proton, and Alpha Monitor (EPAM), the So- 
lar Wind Electron Proton Alpha Monitor (SWEPAM), the Magnetic Field Exper- 
iment (MAG), the Solar Energetic Particle Ionic Charge Analyzer (SEPICA), the 
Solar Isotope Spectrometer (SIS), the Cosmic Ray Isotope Spectrometer (CRIS), 
the Solar Wind Ion Composition Spectrometer (SWICS), the Solar Wind Ions 
Mass Spectrometer (SWIMS), and the Ultra Low Energy Isotope Spectrometer 
(ULEIS). 



2.5.1 SWEPAM 



SWEPAM (McComas et al. , 1998 ) was designed to study the bulk flow and kinetic 
properties of the solar wind, and the suprathermal electrons and ions. Its primary 
goals are to provide solar wind measurements to further our understanding of the 
corona and solar wind formation, as well as the particle acceleration. 

It consists of two main parts, the ion and the electron detector: SWEPAM- 
I and SWEPAM-E, respectively. As ions enter the SWEPAM-I aperture, the 
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Figure 2.11: Schematic of the SWEPAM-E instrument (McComas et al. 1998) 



electrostatic analyser allows only ions of a specific narrow energy and angular 



range to pass through (Figure 2.11). The ions are then detected by channel 
electron multipliers. Similarly to SWEPAM-I, the SWEPAM-E instrument uses 
an electrostatic analyser and channel electron multipliers, but while SWEPAM- 
I separates ions by applying a negative voltage, SWEPAM-E applies a positive 
voltage to filter out electrons of a narrow energy and angle range. The field of view 
is fan-shaped in both instruments with a wide azimuthal angle. The incidence 
measured in the separate electron multipliers determine the angle at which the 
ions or the electrons entered the instrument. The electrons are detected in the 
energy range 1.6 eV - 1.35 keV. The fan shaped field of view covers 95% of Att 
steradian due to the spinning of the spacecraft. 
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Chapter 3 

Coronal hole detection 



Since the discovery of coronal holes, many different methods have been devised to 
improve their detection. From manual drawing to advanced intensity threshold- 
ing techniques, the results can differ considerably. The inconsistency between the 
coronal hole maps is sometimes due to the difference in the detection approaches, 
while in other cases it is an artifact of using different wavelengths that represent 
a different height and temperature range at which the coronal holes are observed. 
In this Chapter we discuss different previous methods developed to detect coro- 
nal holes, and compare the resulting coronal hole maps to those obtained using 
the method developed by the author. The foundation of this thesis research, the 
Coronal Hole Automated Recognition and Monitoring ( CHARM) algorithms, are 
discussed in detail. The diverse use of the detection method is demonstrated by the 
use of several instruments and wavelengths. The research detailed in this Chap- 



ter has been published in the Solar Physics Journal (Krista & Gallagher, 2009). 
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3.1 Introduction 



In order to provide reliable information on their physical properties, CHs must 
be detected in a consistent manner. Initial CH detection methods involved hand- 



drawn maps, such as the method by Harvey & Recely (2002 ) based on the two-day 



average of ground-based He I (10830 A) images and magnetograms. However in 
recent years, a number of automated detection methods have been developed: 



Henney & Harvey (2005) used He I (10830 A) images, magnetograms and the 



above-mentioned hand-drawn CH maps; Malanushenko & Jones (2005) used He I 



imaging spectroscopy to determine CH boundaries; while Toma & Arge (2005) 



have combined space-based EUV images with ground based He I (10830 A) and 
Ha observations and magnetograms to identify coronal hole regions. [Chapman 



& Bromage (2002) and Chapman (2007) constructed synoptic CH maps using 



SOHO Coronal Diagnostic Spectrometer (CDS; Harrison (1997)) measurements. 



Delouille et al. (2007) used a multichannel segmentation method on EUV im- 



ages to differentiate CHs, QS and ARs. More recently, SchoU & Habbal (2008) 



applied an image segmentation technique based on a combination of region and 
edge-based methods using 171 A (Fe ix/x), 195 A (Fe xil), 304 A (He ll) images 
and Michelson Doppler Imager (MDI; Scherrer et al. ( 1995[ )) magnetograms to 
determine CH boundaries. 

The first aspect to be discussed is the visibility of CHs at different wavelengths 
and using different instruments. Some of the above mentioned methods used 
ground-based images which are limited by environmental effects (atmospheric 
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motion can lower the resolution, while clouds disturb the continuity in the ob- 
servations). The analysis and results presented in this thesis are based solely 
on images produced by space-based instruments that are free of environmental 
effects. Some authors used He I images that provide low contrast between CHs 
and the QS. It should be noted though that He I images do have an advantage 
in having a relatively narrow temperature range which can reduce line-of-sight 
obscuration effects. However, EUV and X-ray images have a larger contrast be- 
tween CHs and QS, but are produced with filters of broad temperature response 
that resolve high-temperature solar structures. This can lead to the obscuration 
of CH boundaries, and the effect becomes worse near the solar limb where the 
line-of-sight angle is smaller. 

The second aspect that should be considered is the run-time of the data anal- 
ysis. The algorithm described in this chapter processes a single high-resolution 
1024 X 1024 pixel EUV or X-ray image and a magnetogram under two minutes 
using a computer with a 2.2 GHz Intel Core 2 Duo processor and 2GB 667 MHz 
memory. The detection of CHs in near real-time manner proves a useful applica- 
tion for space weather forecasting purposes. 

The third and most important aspect is the intensity threshold determination 
for CHs. Although determining the intensity threshold is of fundamental im- 
portance to accurately locate CH boundaries (and thus determine their physical 
properties) , many of the methods discussed in the following section rely on arbi- 
trary choices. The method developed by the author determines the CH boundary 
threshold based on an intensity thresholding technique that is also dependent on 
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the mean QS intensity which ensures that the method works rehably at any time 
of the solar cycle. 



3.2 Previous methods of CH detection 

In this section we discuss the methods of CH detection developed by other au- 
thors. As there are numerous approaches, only those detecting CHs on the whole 
visible surface are included here (e.g we exclude studies that only concern polar 
CH identification). The method descriptions concentrate on what instrument (s) 
and wavelength(s) were used, the method of identification, what period was cov- 
ered in the study and whether the CH maps are publicly available. Where it 
is possible, an arbitrarily chosen CH map from the study is compared with a 
corresponding map obtained using the author's method. 



3.2.1 The K. Harvey and Recely method 



Harvey & Recely (2002) used the ground-based National Solar Observatory/Kitt 
Peak (NSO/KP) synoptic rotation magnetograms and He I 10830 A spectrohe- 
liograms to identify coronal holes. The method includes the magnetograms and 
spectroheliograms being remapped into sine-latitude and Carrington longitude, 
but the exact details of the identification process is not discussed. A relevant 



later paper (Henney & Harvey, 2005) indicates that the Harvey & Recely CH 



maps were drawn by hand using computer assistance. 

CHs are known to appear as regions of enhanced emission in He I 10830 A 
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images (Andretta & Jones, 1997). This is due to the overpopulation of the He I 



10830 A transition levels from coronal radiation which leads to increased ab- 
sorption of radiation originating in QS regions, while CHs have reduced coronal 
emission and thus appear as bright regions at 10830 A. To aid the differenti- 



ation of CHs from other QS features Harvey & Recely (2002) defined CHs as 



having less than 20% quiet-Sun network contrast, a predominately unipolar mag- 
netic field containing 75-100% of one polarity, and having the area of at least 2 
supergranules (~109 km^). 

The advantages of the method are the low line-of-sight obscuration in the He I 
10830 A images, and the CHs are observed at their "base" (the chromosphere). 
Also, magnetograms were used to confirm the predominant polarity of the CHs 
detected. The disadvantages of the method are the low intensity contrast between 
the QS and CH regions in the He I images and the hand-drawn quality of the 
CH maps. 



The comparison between the Harvey & Recely (2002) CH map and the corre 



sponding map created by the author is shown in the next section together with 



results of the next method. An extensive database (1989-2002) of the Harvey & 
Recely (2002) CH maps is available onlin^ 



3.2.2 The Henney and J. Harvey method 



Henney & Harvey (2005) developed an automated CH detection algorithm us- 



ing the above mentioned hand-drawn maps, SOLIS Vector SpectroMagnetograph 



ftp: / /nsokp. nso.edu /kpvt / synoptic / choles / 
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(VSM) helium spectroheliograms and photospheric magnetograms. As a first 
step, consecutive pairs of spectroheliograms were derotated to the time of the 
later observation and remapped into sine-latitude and longitude. This was fol- 
lowed by taking two-day averages of the image pairs, after which the previously 
mentioned hand-drawn CH contours were overlaid on the average magnetogram 
and helium spectroheliogram in order to determine the normalized He I 10830 A 
equivalent width and mean percent unipolarity. 

The next step involved taking only values above 1/10 of the median of the 
positive values in the average helium image. This was followed by removing 
the negative values, smoothing the image, and all values were set to a fixed 
number. Next, a square kernel function was used as the shape-operator to fill 
gaps and connect nearby regions while small regions were removed (the same 2 



supergranular size threshold was applied as in Harvey & Recely (2002)). The 
produced image was then used as a mask to set coinciding pixels in the median- 
selected image to a high value. The final image was then obtained by smoothing, 
and removing small CHs and those of low unipolarity percentage. 

The advantages of the method are the automated quality of the detection, 
the use of a wavelength that provides low line-of-sight obscuration, and using 
magnetograms to confirm the predominant polarity of the CHs detected. The 
disadvantages are the low contrast between the CH and QS regions that could 
impair the detection of the CH boundary, and the intensity threshold, which was 
based on the hand-drawn maps. 

It must also be noted that the publicly available SOLIS CH maps are Car- 
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rington rotation maps which are created from stacked central meridian plots over 
~27 days. "Time-delayed" Carrington plots are only an approximation of the 
global CH locations and morphology. (I.e. by the time the Carrington map is 
created one end of the image is almost a month "older" than the other end of the 
image.) 

The [Henney &: Harvey (2005) Carrington rotation CH maps are available from 
2004, and are continuously updated. The database is available onlin^ 



Figure 3.1 shows a NSO/KP He 1 10830 A spectroheliogram and a SOHO/EIT 



195 A image taken on 14 July 2003 (top left and right, respectively), and the com- 



puter assisted hand-drawn maps by Harvey & Recely ( 2002 ) and the correspond 



ing multi-step automated CH detection by Henney & Harvey (2005) (bottom left 
and right, respectively). For comparison, the CH contours obtained using the 



authors' method is shown in Figure 3.2 on a Lambert equal area projection map. 



3.2.3 The Chapman and Bromage method 



Chapman & Bromage ( |2002 ) used synodic data from the SOHO/CDS Normal 
Incidence Spectrometer (NIS) to develop a partially automated CH detection 
method. Nine 4x4 arcmin rasters were used to create an image-strip along the 
solar meridian. The strips were obtained daily to produce a Carrington map. 
Among the available NIS wavelengths the Mg ix line observations were used for 
identification, since in these images CHs are distinguishable as dark regions. After 
the low emission regions were extracted, the CHs were manually differentiated 



^ http://solis.nso.edu/vsm/vsm_maps.php "Archive of solar wind source synoptic maps" 
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Figure 3.1: The top left and right images are a NSO/KP He I 10830 A spectro- 
heliogram and a SOHO/EIT 195 A image taken on 14 July 2003, respectively. The 
bottom left and right images are the CH maps created with the computer assisted 



courtesy of (Henney Sz Harvey, 2005). 



hand-drawn method devised by Harvey & Recely (2002) and the corresponding 



map from the Henney & Harvey (2005) automated CH detection method. Figure 



from other low intensity regions using He I 10830 A images. Carrington rotation 
CH maps were then generated by joining raster images from consecutive days, 
removing overlapping areas. The method was applied to solar images in the 



period July 1997 - Sep 2001. Chapman & Bromage (2002) found that the CH 



area cycle was anti-correlated with the sunspot cycle. 
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EIT Fe XII (195 A) 13-Jul-2003 23:24:10,899 
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Figure 3.2: The SOHO/EIT CH map created by the author to compare with 
the results of the [Harvey &: Recely and Henney &: Harvey methods as seen in 
r igure O An EIT 195 A observation taken on 13 July 2003 and a corresponding 
magnetogram was used in the detection. (The numbers represent different groups 
of CHs, as discussed later in the chapter.) 



The method was further developed in the following years. In the thesis re- 



search of Chapman (2007), Mg x/Mg ix temperature ratio maps were used in 



conjunction with Mg X and Mg IX intensity images. Upper and lower boundaries 
were set to specify which low temperature and low intensity regions classified 



as CHs. An example CH map by Chapman showing the 13 July 2003 large 



transequatorial CH (at Carrington longitude 180-90) is shown in Figure 3.3 (for 



comparison see Figures 3.1 and 3.2 
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Figure 3.3: The CH map created by Chapman (2007) for Carrington rotation 
2005. 



The advantages of the method are the automated quahty of the detection, the 
abihty to investigate the temperature of the low intensity regions and the narrow 
bandpass imaging which hmits hne-of-sight obscuration effects. The shortcom- 
ings are the time-consuming analysis, the time-delayed raster-stacked Carrington 
maps and the absence of magnetic field analysis within the low temperature and 
intensity regions. 

The analysis was completed for the period 1998-2007, and the CH Carrington 
rotation maps are available onlin^ 

3.2.4 The Scholl and Habbal method 



Scholl fc Habbal (2008) used SOHO/EIT 304 A, 171 A, 195 A and 284 A observa- 



•'http://www.star.uclan.ac.uk/~sac/coronalJioles/images.html 
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tions and MDI magnetograms to identify CHs. A combination of image processing 
techniques were used to differentiate CHs and filaments. Although the steps are 
not detailed in the corresponding article, a list of image processing operations 
were given. Due to their lower contrast, the 304 A images were processed in a 
different way from the EUV images. From the numerous image processing steps, 
only the main ones are listed here. After calibration and normalisation, limb- 
darkening was removed from the 304 A images, and the contrast was enhanced 
in the images of all four wavelengths. This was followed by smoothing to connect 
dark regions, image-segmentation based on histogram-equalisation, and region- 
growing to isolate dark regions. A contouring method was used to extract regions 
of interest and remove small features. The candidate regions were then obtained 
separately at the four wavelengths. Finally, the regions were classified as CHs or 
filaments based on their magnetic field profile and the threshold values specified 
by the authors. 



An example for the SchoU & Habbal (2008) method is shown in Figure 3.4 



with CHs (orange, right panels) and filaments contoured (black, middle panels). 
Filaments (green) and ARs (red) from the Solar Feature Catalogue (SFC) are 



shown for comparison (Bentley et al. 2004). The top, middle and bottom row 
maps are based on 304 A, 171 A, and 195 A SOHO/EIT observations of 27 May 
1997. In the right panels the NSO/KP CHs are overplotted in black contours for 
comparison. The corresponding CH map determined by our algorithm is shown 



in Figure 3.5 



The advantage of this method is the automated approach, the use of several 
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Figure 3.4: The feature classification map by 



Scholl k Habbal 



(2008) showing 



CHs (orange, right panels), filaments (black, middle panels). SFC filaments (green) 
and active regions (red) are shown for comparison (Bentley et al.. 2004). The top. 



middle and bottom row maps are based on 304 A, 171 A, and 195 A SOHO/EIT 
observations, respectively. Observations were made on 27 May 1997. NSO/KP 
coronal holes are overplotted as black contours for comparison. The blue line 
separates the opposite-polarity magnetic-field areas shown as gray and white in 
the left and middle panels. 



wavelengths which map CHs from different perspectives of temperature and atmo- 
spheric height. Also, magnetograms were used to differentiate between CHs and 
filaments. The disadvantages are the time-consuming analysis of the numerous 
images, and the inconsistent success rate of the method during the solar cycle 
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EIT Fe XII (195 A) 23-Mar-2003 01:13:42.909 




Longitude (degrees) 

Figure 3.5: The SOHO/EIT CH map created by the author to compare with the 
results of the lScholl fc Habball An EIT 195 A observation taken on 23 March 2003 
and a corresponding magnetogram was used in the detection. 

(it was shown to be more successful during the solar minimum). It is unclear 
whether an online database is available and what period was processed. 



3.3 The Ci/^i?M algorithm 



The CHARM algorithm presented in this section is a fast and automated method 
which allows the user to determine CH boundary thresholds in a consistent man- 
ner throughout the solar cycle. The detection tool can also be used at different 
wavelengths and with different instruments. Additional packages allow the user 
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to study the short-term changes in CH boundaries, and to determine the ar- 
rival time of high-speed solar wind streams at Earth based on the CH boundary 
locations. These are discussed in Chapter |4] and Chapter [5| respectively. 

3.3.1 Data preparation 

The detection tool requires the user to specify the wavelength and instrument. 
The use of EUV 195 A and X-ray wavelengths are advised for CH detection as 
they provide better CH boundary contrast than e.g. 171 A, 284 A 304 A or He I 
10830 A images. The start time and duration of the analysis period are also de- 



termined by the user. This is followed by standard SolarSoft procedures Freeland 



& Handy (1998) to calibrate the EUV and X-ray images. These include: 

1. Dark current subtraction. Taking a time exposure with the shutter closed 
produces the thermal signature known as "dark count", caused by thermal elec- 
tron excitation due to the temperature of the device. This is the zero flux CCD 
response, which is subtracted from the raw image. 

2. Instrument grid subtraction. The aluminum filter near the focal plane of the 
instrument casts a grid-like shadow on the CCD detector, which has to be sub- 
tracted from the image. 

3. CCD response correction, also known as "flat flelding" . Each CCD pixel has a 
different sensitivity to incoming photons. This can be corrected for by dividing 
the data by flat-fleld exposures taken before the observation. 

4. Exposure time normalisation ensures that each observation is normalised ac- 
cording to the exposure time of the observation. The resulting image intensity 
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unit is data numbers per second (DNs; the measure of photon incidence in the 

CCD camera). 

5. Filter normalisation corrects for the variable transmittivity of the clear and 
aluminium filters. 

The output is an index structure containing the observation details (stored in 
the header) and an image data cube. In this thesis work the calibration routines 
"eit_prep", "secchi.prep" and "xrt_prep" were used to prepare level-1 data. The 
MDI magnetograms and the PLASTIC and ACE solar wind data needed no 
further calibration. 

3.3.2 Data selection 

The EUV images are used to detect low intensity regions, which are complemented 
by magnetograms to determine whether a region is of a predominant polarity 
(typical of CHs). SOHO provides high resolution (1024x1024 pixel) full-disk 
EUV images and magnetograms, and has the largest database including over 10 
years of data. For this reason, the present thesis research is mainly based on data 
from the SOHO EIT and MDI instruments. Hinode/XRT and STEREO/EUVI 
data are used for comparison and confirmation. 

The magnetograms analysed can be chosen to be one-minute or five-minute 
averaged full-disk magnetograms. In this work we used the five-minute images, 
as these are the lowest noise images produced daily. The per pixel noise in 
the Level 1.8 calibrated magnetograms is ~30 G for a 1-minute magnetogram 
and ~15 G for a 5-minute magnetogram. The Level 1.8 magnetograms are the 
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definitively calibrated MDI images of the solar line-of-sight magnetic field. 

In order to detect CHs at a given time, a magnetogram is found closest to 
the required time. Next, the closest EUV image is selected within ±30 minutes 
(this time range can of course be changed). The reason for selecting the MDI 
magnetogram before the EUV image is the smaller number of MDI observations 
per day compared to EIT 195 A observations. There are fifteen MDI magne- 
tograms created in a day and often over a hundred EUV images. Next, the 
intensity threshold is determined for the low intensity regions (as explained in 



Section 3.3.4), which is followed by the magnetogram being rigid-rotated to the 
time of the EUV observation to match the corresponding locations on the solar 
disk. The magnetograms are then used to determine whether the low intensity 
regions are CHs. 

3.3.3 The Lambert cylindrical equal-area projection 

The Lambert cylindrical equal-area projection allows a spherical surface to be 
projected onto a cylindrical surface whilst preserving the area and producing 
parallel latitudinal (and longitudinal) lines within approximately -60 and +60 



degrees latitude and longitude McAteer et al. (2005). This projection was chosen 



to aid the local intensity thresholding procedure, which uses "window-stepping" 



across the image (Section 3.3.4). It also provides a straightforward determination 
of the CH area and the boundary positions. 

The cylindrical equal-area projection was the fourth of seven projections in- 
vented by Johann Heinrich Lambert in 1772. The equations that describe the 
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Figure 3.6: The Tissot indicatrix showing the longitudinal and latitudinal line 
positions and the shape distortions after the Lambert cylindrical equal-area pro- 
jection of the world map. Figure courtesy of Eric Gaba and U.S. NGDC World 
Coast Line. 



projection are 



X = R{X — Ao) cos (ps 



(3.1) 



y = Rsmcpl cos( 



(3.2) 



where R is the sphere radius, Aq is the central meridian location (set to 0° in 
our projection), A is the longitude, is the latitude, and (ps is the standard 



parallel (set to 0° in the equator centered projection), as described by Snyder 



(1987). As an example, the Lambert projection of the whole Earth's surface is 



shown in Figure |3.6[ Note that the longitudinal lines are equally spaced, while 
the latitudinal lines are sine spaced. 

The algorithm used for the projection is lambert .project (a modified version 
of map .carringtonS) . The Lambert cylindrical equal- area projection is better 
known as Carrington projection in heliographic work. However, in the majority 
of cases Carrington maps are stacked maps of solar central meridian daily strip 
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images. The Lambert projection procedure projects the visible solar surface into 
a cylindrical map, however, uncertainties are introduced as the surface visibility 
drops near the solar limb. For this reason the projection is not as reliable above 



±60° latitude and longitude. Figure 3.7 shows the projection deformation in 



latitudinal and longitudinal lines, which are drawn at 15° intervals. By default, 
the algorithm limits the projection to ~ ±80° latitude and longitude due to lack 
of information near the limb. In addition to this, ofF-limb regions are visible in the 
Lambert map corners with varying scale due to the yearly variation is the solar 
Bo angle. This is caused by the inclination of the Earth's orbit to the ecliptic 
giving a larger view of the northern or southern solar pole during an orbital year 
of Earth. This variation is corrected for in the algorithm. As a result, the solar 
equator is located in the image center in all Lambert maps. The off-limb corners 
are excluded from all analysis. 

3.3.4 Detecting low intensity regions 

Intensity thresholding is a powerful method to distinguish regions of different 
intensities in an image (e.g. CHs and QS regions). The intensity histogram 
of an image gives a multimodal distribution where each frequency distribution 



corresponds to a feature in the image Belkasim et al. (2003); Gallagher et al. 



(1998 1999a); Gonzalez & Woods (2002). In the case of 195 A images, each 



major feature (ARs, low intensity regions - LIRs, and QS) corresponds to a 
distribution in the intensity histogram. When an image is dominated by QS, 
the corresponding intensity histogram shows a unimodal distribution (top panel 



102 



3.3 The CHARM algorithm 



EIT 1-Jun-2008 01:14, 15° grid, B0 = -0.624135 




X [arcsecs] 



EIT 1-Jun-2008 01:14, 15° grid 



- \^ 


1^ 


■ 




1 


■ 


r r ri ~ 
















LZ 




















/ 


1,1 : 






















J 












L. 


I 


I \ 




1 

Ml 


m 






m 


■ 


■LA 



-75 -60 -45 -50 -15 1 5 50 45 60 75 
Longitude (degrees) 

Figure 3.7: The EIT disk image (top) and tlie Lambert cylindrical equal area 
projection of the visible solar surface with 15° grid spacing in latitude and longitude 
(bottom). (Observation: SOHO/EIT 195 A, 1 Jan 2008). 



103 



3. CORONAL HOLE DETECTION 



SECCHI A 2-Jan-2008 00:05 

90 I ' ~ 



ntensity histogram 



a; 30 h 

3 

dJ 

3 30 
-90 



-45 45 

Longitude{ degrees) 

SECCHI A 15-Jan-2008 00:05 
90 




-45 45 
Longitude{ degrees) 




30 40 60 80 100 
I(DNS) 



Intensity histogram 




20 40 60 80 100 

I(DNS) 



Figure 3.8: Top: EUVI 195 A Lambert map showing no CHs, and the cor- 
responding intensity histogram with a unimodal distribution. Bottom: Lambert 
map showing a large CH, and the corresponding intensity histogram with a bi- 
modal distribution. 



of Figure 3.8). When a relatively large part of the image is covered by LIRs, 



the resulting intensity histogram shows a bimodal distribution (bottom panel of 



Figure 3.8), where the local minimum corresponds to the intensity threshold of 



the LIR boundary. Figure |3.8| gives the intensity measured in DNs as a function 
of CCD pixel number frequency. 

The minimum between the LIR and QS distribution was enhanced using a par- 
titioning operation. By partitioning the solar image to sub-images, we obtained 
local histograms with more defined minima which aided the determination of 
thresholds. The sub-images used were 1/2, 1/3, 1/4 and 1/5 of the original Lam- 
bert map size. Each of the sub-image sizes were used to step across the map and 



obtain local intensity histograms from overlapping sub-images (Figure 3.9). This 
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Figure 3.9: A schematic of the Lambert map being partitioned into overlapping 
sub- images. 

process provided more clearly discernible local minima for thresholding. This 
can be explained by different size sub-images having different proportions of QS 
and LIR and hence frequency distributions. The threshold minima are defined as 
local minima that are located within 30-70% of the QS intensity and are wider 
than 6 DNs. The threshold range was set based on CH boundary intensity tests 
through the solar cycle. The range is needed to reduce run time and to exclude 
occasional low intensity features inside CHs, or high intensity features in the QS. 

As mentioned above, the intensity range where the thresholds were identified 
depends on the QS intensity. For this reason the method is reliable at any time of 
the solar cycle regardless of the change in the overall intensity. It should be noted 
however, that with increasing solar activity, the acquired CH boundaries are less 
accurate due to bright coronal loops intercepting the line-of-sight and obscuring 
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parts of the CH boundary. 

Once the intensity thresholds are obtained from the sub-images the final in- 
tensity threshold is determined by taking the most frequent value. The threshold 
can then be used to contour LIRs on the same Lambert map. In our studies we 
chose to determine the threshold for five images in a day and set the median value 
of the thresholds as the final threshold. This is done to overcome small intensity 
variations in time and secure a more reliable threshold. An example is shown 



in Figure 3.10 for a CH complex observed with four different instruments in two 
wavebands (195 A and X-ray). The observation times are different to ensure that 
the CH is seen at the same central location observing with different satellites. 



The threshold histograms (right column of Figure 3.10) show distinct threshold 
peaks, which are used to contour LIRs. We would like to draw attention to the 
large-scale similarity in CH morphology despite the use of different instruments 
and observation times. This similarity is aided by the CH being close to the 
disk center, where line-of-sight obscuration is minimal. However, it is important 
to note that in other cases differences may appear between EUV and X-ray CH 
observations due to the difference in the temperature range they observe at. Ob- 
servations made at different angles from the Sun (e.g. SOHO, STEREO A and 
B are at different orbital locations) may also lead to differences in the CH area 
and boundary location due possible line-of-sight obscuration. 
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Figure 3.10: Left STEREO A/EUVI, SOHO/EIT, Hinode/XRT and STEREO 
B/EUVI Lambert maps with LIR contours. The images show a complex CH struc- 
ture approximately at the same location on the map. Right: Threshold histogram, 
where the peak value is the optimal threshold. 
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3.3.5 Differentiating coronal holes from filaments 

CHs were distinguished from filaments using MDI magnetograms. As previously 
mentioned, the magnetograms were chosen first in the analysis based on the 
time of interest. The magnetogram was then differentially rotated to the EUV 
image observation time. The LIR contours obtained from the EUV image was 
then overplotted on the magnetogram and the flux imbalance was determined 
for the contoured regions. CHs are well known to have a dominant polarity, 



while filaments have a bi-polar distribution of polarities SchoU & Habbal (2008). 



Figure 3.11 shows examples of LIRs observed on 15 June 2008 and 22 March 
2007, respectively. The 15 June 2008 LIR is an example of a CH with a neg- 
ative dominant polarity, while the 22 March 2007 LIR is a CH with a positive 



dominant polarity (first and second row in Figure 3.11). The skewness of the 
field strength distribution defines the scale of the flux imbalance in the observed 
region of interest. For the two CHs the skewness was found to be considerable 
with values —2.9 and 6.4, respectively. The LIR highlighted on the bottom panel 



magnetogram in Figure 3.11 was found to be a filament based on the very small 



value of skewness (0.2) in its line-of-sight magnetic field strength distribution. 
3.3.6 Coronal hole grouping 

The grouping of CHs was done mainly for high-speed forecasting purposes (Chap- 
ter |5]). However, it also simplifies the tracking and referencing of complex CH 
groups. The group association is not founded on physical concepts; it is simply 
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Figure 3.11: Left: EIT Lambert maps with LIR contours. Middle: MDI Lam- 
bert maps witli a LIR contours overlaid. Right: Blqs intensity histogram of the 
contoured regions. Regions highlighted in the first and second middle panels are 
CHs {^LOS histograms skewed), the third middle panel shows a filament {Blqs 
histogram not skewed). 

a measure of distance between CHs. 

Each CH is treated as a polygon located on the cylindrical Lambert map. 
The pixel coordinates of the polygons are converted into heliographic coordinates 
after which the Euclidean distance is determined between all coordinate points. 
The smallest distance between the two polygons' coordinate points is the mini- 
mum distance which determines whether two CHs belong to the same group. The 
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Figure 3.12: Top: The LIRs numbered are CHs, the white contours with no 
numbers are LIRs that did not meet the requirements to quahfy as CHs (i.e. due 
to too smah area or not large enough flux imbalance). Bottom: CHs after the 
grouping procedure. CHs of the same group have the same number and colour. 

distance threshold is set to 10°, above which CHs are not grouped. 



Let us consider the LIRs detected in the EUV image shown in the Figure 3.12 
top panel. These are numbered from 1 to 12 in the descending order of their areas 
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Figure 3.13: The CHs selected from the LIRs. 

(i.e. the largest CH is numbered 1). During the magnetic field analysis nine of 
the LIRs, numbered 1, 2, 4, 5, 6, 8, 10, 11, 12, were classified as CHs. The list of 



these CHs is shown in Figure 3.13 These CH numbers are non-consecutive due 
to the discarded LIRs. 

First, the CHs are grouped with each other based on their distances as de- 



scribed before. The first and second column in Figure 3.14 are the CH numbers 
in order, and the neighbouring CHs found at less than 10° distance, respectively. 
Repetitions, such as CH 4 paired with CH 5, and CH 5 paired with CH 4, have 
to be removed from the list. Note that the repeats with the higher CH number 
(thus the smaller area) are the ones to be removed. This ensures that the CH 
with the largest area becomes the Group Leader, and the smaller neighbouring 
CHs become the Group Members. 

The problem becomes complicated when a CH (e.g. CH 11) is not directly 
associated with a Group Leader (i.e. CH 4), but with a Group Member (e.g. 
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Figure 3.14: The CHs and their associated pairs. Each CH has beed given a 
different colour to highhght the pairs. The repetitions (crossed out with a red hne) 
are removed from the hst. 



CH 5) which is associated with a Group Leader (CH 4). Because all three CHs 
( CH 4, 5 and 11) are close to each other, they become grouped, and the largest 
CH's number will be the Group Number (4 in this case). 



Figure 3.15) shows the previous table, but with repetitions excluded and fur- 
ther associations included (e.g. CH 11 is a member of Group 4 through being 
close to CH 5). Next, the Group Leaders are given consecutive group numbers 



{Final Groups), as shown in the table in Figure 3.16 (also see the corresponding 



CHs in Figure 3.12 bottom panel). 
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Figure 3.15: The CHs and their associated pairs. 
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Figure 3.16: The CHs and their associated pairs were renumbered to produce 
consecutive Final Group numbers. 



3.4 Conclusions 

This Chapter introduced the Ci/^i?M algorithm, which was developed to identify 
CHs in fast and automated manner. It allows the use of different instruments 
and wavelengths (X-ray and EUV), and gives consistent CH boundary thresholds 
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throughout the solar cycle. 

Previous methods have used arbitrary thresholds or a computationally inten- 
sive approach, while the CHARM method only uses a single EUV image and a 
magnetogram to detect CHs at the time of interest. The segmentation algorithm 
first identifies LIRs using local intensity histograms of either 195 A or X-ray 
observations, and then exploits MDI magnetograms to differentiate CHs from 
filaments by the flux imbalance determined for the identified regions. This was 
achieved assuming that CHs have a predominant polarity as opposed to filaments. 
Once the CH regions were determined, the position, area and mean -Blos can be 
calculated for multiple CHs. Finally, the CHs are grouped mainly for forecasting 
purposes. However, this can also aid the tracking and referencing of CHs. 

The main limitation of the method results from near-limb effects, which in- 
clude: line-of-sight obscuration in EUV and X-ray images, increase of noise in 
-Blos measurements, and interpolation effects caused by the Lambert equal-area 
projection. These near-limb effects can cause incorrect boundary determination, 
the inability to determine a dominant polarity, and incorrect area determination. 
Also, some CH boundaries located near the limb were found to have different 
threshold intensities which depended on the degree to which they were obscured. 
For this reason, only CHs within -60 and +60 degrees latitude and longitude were 
included in this analysis. 

Due to the fast and automated approach of the present method, it allows the 
analysis of large databases (e.g. SOHO database). This provides an opportunity 
to study the long-term evolution of CHs (Chapter^. On a smaller scale, the 
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method can also be applied to analyse the small scale changes in CH boundaries, 
which may provide information on the processes governing CH growth and decay 
(Chapter g. 
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Chapter 4 

Short-term evolution of coronal holes 



The interaction of open and closed field lines at CH boundaries is widely accepted 
to be due to interchange magnetic reconnection. To date, it is unclear how CH 
boundaries vary on short timescales and what velocities this occurs at. The author 
developed an automated boundary tracking method to study short-term boundary 
displacements, which were found to be isotropic and have typical speeds of 0- 
2 km s^^. The magnetic reconnection rates at CH boundaries were calculated to 
be 0-3x10^^, and using the interchange reconnection model the diffusion coef- 
ficients were found to be 0-3x10^^ err? s~^ . These results are consistent with 
interchange reconnection in the corona, driven by the random granular motion of 
open and closed fields in the photosphere. This work has been accepted for pub- 



lication in Astrophysical Journal Letters (Krista, Gallagher & Bloomfield, 2011). 
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4.1 Short-term changes in CH boundaries 

Small-scale changes in the CH morphology are thought to result from magnetic 
reconnection occurring between CH open fields and small closed fields found in the 



QS as well as in CHs (Madjarska et al. , 2004). Using the Coronal Hole Evolution 



CHE VOL algorithm we investigate the driving force behind the open and closed 
field interaction and determine the magnetic reconnection rate and the diffusion 
coefficient based on the observed displacements in CH boundaries. 

A number of authors have previously studied the large-scale and long-term 
evolution of CHs and established the changes in the total CH area, flux, the 
latitudinal distribution over the solar cycle and the key role CHs play in the 



dipole-switch of the Sun (Bravo &; Stewart 1997 Chapman & Bromage, 2002 



2002 


Wang 


2009 


Wang et al. 


1996) 



scale changes in CH boundaries have only been studied in detail in recent years 



using high- resolution observations (Madjarska & Wiegelmann 2009 Subrama- 



nian et a/., 2010), it was already suggested in the late 1970s that the small-scale 



and short-term evolution of CH boundaries is strongly connected to bright point 



occurrences (Davis, 1985 Nolte et al, 1978) and that reconnection could be a 



means of flux-transportation and "diffusion" of the open flux (Nolte et al, 1978). 



A more in-depth discussion on reconnection occurring at CH boundaries was pro- 



vided by MuUan & Ahmad (1982), who proposed that reconnection was also a 
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driver of the solar wind emanating from CHs. 



4.2 Interchange reconnection 

The process of continuous reconnection between open and closed field lines is 



2010; 


Fisk, 


2005 


Raju 



et al. 2005 Wang & Sheeley, 1993 2004 Wang et al. 1996). As suggested by 



Fisk & Schwadron (2001) and Fisk (2005), interchange reconnection allows the 



transfer of magnetic flux from one point to another, which can be considered 
as a diffusion-convection process resulting from convective motions and random 
motions caused by non-stationary supergranular flows. The equation for this 
diffusion process is based on the induction equation which we derive using the 
Maxwell equations: 



^ X. . 1 



V X E 



V-E 



P 



So 

V ■ B = 



(4.1) 
(4.2) 
(4.3) 
(4.4) 



where B is the magnetic field strength, E is the electric field strength, is the 
magnetic permeability, j is the current density, c is the speed of light, p is the 
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charge density and Eq is the electric permittivity. We also include the Ohm 



equation: 



j = a(E + vxB) , 



(4.5) 



where a is the electrical conductivity, and v is the plasma velocity. In the MHD 
approximation the plasma velocity is much smaller than the speed of light {v << 



c), thus VxB>>l/c dEi/dt, and hence Equation 4.1 reduces to 



V X B = /ij 



(4.6) 



which we rewrite as j = l//i (V x B) and substitute into Equation 4.5 to get 



E = -V X B + —(V X B) 
/icr 



(4.7) 



This equation is substituted into Equation 4.2 with rj = 1/fia to get 



dB 

~dt 



V X (v X B) - V X (r^V x B) 



(4. 



which reduces to 



dB 

~dt 



V X (v X B) - r^V^B 



(4.9) 



also known as the induction equation, where the first term on the right is the 



convection term, and the second is the diffusion term. Here, we assume that the 
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magnetic field lines diffuse due to random motions in the medium. The resistive 
diffusion is ignored (V x (v x B) = 0) as we only consider the field line motions 
due to the velocity of the medium. If only the radial motions along the solar 



surface are considered, Equation 4.9 can be written as 



dB 

^ = -V, ■ (Vrand^,) , (4.10) 



where Vrand is the plasma velocity along the surface due to random motions, and 
Bz is the radial component of the magnetic field. Following Fisk & Schwadron 



(2001) it can be shown that 



^ = -WI{kBz) , (4.11) 



where B^ is the mean radial magnetic field and k is the diffusion coefficient tensor, 
which in a simplified form can be written as 



where 5r is the open magnetic field displacement distance along the surface and 
5t is the time elapsed. 

Despite all of the flux transfer processes involved being convective in nature, 
the non-stationary random supergranular convection is included in the diffusion 
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term while those that are uniform (such as differential rotation and meridional 
flow) are included in the convection term. As supergranular motions are too slow 
for flux transportation, interchange reconnection was suggested as a faster mech- 
anism for field-line diffusion. Interchange reconnection allows for the transfer of 
magnetic flux from one point to another through magnetic reconnection between 



open and closed field lines, which Fisk & Schwadron (2001) modelled using a 
diffusion-convection equation, 

= kVIB, - V ■ (v^onvB,) , (4.13) 

where Vconv is the uniform convective flow velocity. It is assumed that the radial 
component of the magnetic field is equivalent of a surface density and is subject 
to diffusion and convection along the solar surface. The diffusion term corre- 
sponds to random supergranular motions, while all convective motions including 
the differential rotation and meridional flows are included in Vconv Because su- 
pergranular motions are considered too slow to transport open magnetic flux on 
the Sun. Instead, a faster process of interchange reconnections between numerous 
open and closed field lines is considered as a process for diffusion. 



A graphical representation of interchange reconnection is given in Figure 4.1 
where a QS loop has moved close to a CH boundary due to random granular 
motions. Magnetic reconnection between the open and closed field leads to the 
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Figure 4.1: Interchange reconnection: Tiie location of the CH and the QS are 
indicated on the cartoon with a hne representing the boundary between the two 
regions. The scenario on the left shows a QS loop located near a CH open field 
which yields magnetic reconnection due to the proximity of the opposite polarity 
loop footpoint to the open field. The result of the reconnection is shown on the 
right - a small loop is created inside the CH with opposite magnetic orientation 
and the open field line "relocated" from within the CH to the further foot-point of 
the previously larger loop. This process can explain the rigid rotation as well as 



the growth of CHs (modified from Fisk, 2005) 



displacement of the open field and thus the CH boundary moves further out. 
During this process the larger loop is replaced by a smaller loop with oppositely 
oriented footpoint polarities imbedded inside the CH, while the open field line is 



"relocated" to the former foot-point of the larger loop. 



Interchange reconnection at CH boundaries is also supported by observational 
results, e.g., the prevalence of very small loops inside CHs with heights <3 Mm 



and larger loops with heights >3 Mm outside CHs (Wiegelmann & Solanki, 2004), 



the presence of bidirectional jets (Doyle et a/. 2006b Kamio et a/. 2009 Mad- 



jarska et ai, 2004), and the boundary displacements observed due to the emer- 
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gence and disappearance of bright small-scale loops in the form of bright points 



(Madjarska & Wiegelmann 2009; Subramanian et al. 2010). 



In this Chapter, we study the short-term evolution of CH boundaries by ob- 
jectively determining CH boundaries based on intensity thresholding methods, 
as well as introducing a new automated boundary tracking tool (CHEVOL) that 



allows us to follow the changes that occur (detailed in Section 4.3). In Section 



4.4 the magnetic reconnection rate and diffusion coefficient are determined using 



the interchange reconnection model, and the measured boundary displacement 
velocities. The processes that influence the CH boundary evolution are discussed 
in Section 1431 



4.3 Observations and Data Analysis 

Finding objective CH boundary intensity thresholds is fundamental in deriving 
consistent physical properties (e.g., area, magnetic field strength) throughout the 
CH evolution. For this reason we used the CHARM algorithm (see Chapter [sj 
on a localised small region of interest. The size of the studied region was chosen 
to be a 31° x 35° "box", which was centered on each analysed CH. 

In this study we used data most suitable for CH detection: 12 minute cadence 
195 A images from SOHO/EIT to locate low intensity regions, and 96 minute 
cadence magnetograms from SOHO/MDI to detect predominant polarities and 
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hence identify CHs. The magnetograms were also used to determine the magnetic 
properties of the CHs and the surrounding QS. The study includes eight relatively 
small CHs, which were chosen arbitrarily between 2008 September 8 and 2009 
October 17, and tracked for two to six days. The CHs studied included both 
predominantly positive and negative polarity holes, with equal numbers located 
in the northern and southern hemispheres. All of the CHs studied were non- 
fragmenting and of relatively simple morphology throughout their evolution. 

The CHE VOL method only requires the start-date and duration of the anal- 
ysis, and the location of the region of interest to be given. The start-date should 
preferably be given for when the studied CH is in the eastern hemisphere. This 
allows the CH to be studied for a longer period of time before it rotates off-limb. 
To avoid near-limb effects, the CH should not be further than 50° from the central 
meridian. The location of the CH is provided by the user by clicking on the EIT 
Lambert map that appears after the program is executed. The user is requested 
to click south-east of the studied CH, as this determines the lower-left corner of 
the box. The box location is saved together with the start-date, thus the region 
specification only needs to be provided once for each CH studied. The size of 
the box can be altered if needed. Prom this point on the analysis is completely 
automated, and the specified box "follows" the CH in consecutive EIT images in 
the time period requested. 

All EIT images are initially transformed into Lambert equal-area projection 
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Figure 4.2: CH contours in an EIT 195 A Lambert projection map (a) and the 
corresponding sub-image (b). MDI sub-image with the CH contours (c) (>50 G 
fields: yeUow, <-50 G fields: blue) and the corresponding -Blos histogram (d) 
(magenta: CH, black: QS). 



maps to aid the intensity thresholding and to simphfy the tracking of the region 



of interest. Figure 4.2 1 shows the EIT Lambert projection map with the CH 



contours and the box around one of the CHs studied (only the largest CH in 



the box was a candidate in our study). Figures 4.2 d and 4.2: show the EIT and 



MDI close-up images of the box with the CH contours, while Figure 4.21 shows 
the histogram of the line-of-sight magnetic field (-Blos) from the largest CH (ma- 
genta line) and the surrounding QS (black line), excluding small satellite CHs. 
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The skewness of the -Blos histogram is used to estabhsh whether a predominant 
polarity is found, thus allowing the differentiation of CHs (skewness > 0.5) from 
other low intensity regions (skewness < 0.5). The agreement between the visual 
and the analytically-obtained CH boundary, together with the predominant polar- 
ity found in the contoured regions, validate the appropriateness of the boundary 
positions and justifies the CH status of the detected regions. In principle, the 
intensity of a CH boundary can change due to the line-of-sight obscuration at 
higher longitudes and latitudes. However, all CH candidates in this study lie 
within ±50° in longitude and latitude, which lowers both the obscuration effect 
in EUV images and the errors in line-of-sight magnetic field measurements. 
Once a CH was identified, consecutive boundary observations were compared. 



Considering an example of contour pairs as shown in Figure 4.3 1, the centroid of 
the earlier observed contour (solid line, observed on 2008 September 8 at 05:48 
UT) was considered as the origin of the frame of reference for both contours. The 
later contour (dashed line, observed on 2008 September 8 at 17:00 UT) will have 
moved with the synodic solar rotation rate compared to the earlier contour, due 
to solar rotation. For this reason a correction was applied to the second contour 
coordinates. The locations of both contours were measured radially from the 
earlier centroid in 5° segments and anti-clockwise from west of the earlier con- 
tour's centroid (Figure |43b). Each boundary segment location was determined 



as the median of the contour point distances from the reference centroid. The 
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displacement of a segment {5r) and the time difference between the two obser- 
vations yielded the displacement velocity of the segment. This calculation was 
repeated to get all the boundary displacement velocities over 360° (Figure [43):) ■ 



(Note that the velocity mentioned here depends of the temporal resolution of the 
used instrument and the time-scale at which the boundary displacement hap- 
pens.) The displacement velocities measured for the whole observation period of 



a single CH were also determined, as shown in the 2D histogram of Figure 4.3 i. 
Here, dark regions correspond to more frequent occurrences of the boundary ve- 
locity values, where the observations were binned in 20° and 0.4 km s~^ intervals. 
These calculations were further extended to include every boundary displacement 
velocity of each CH in the study, as detailed in the following Section. 

4.4 Results 

4.4.1 Coronal hole boundary displacement velocities 

In the case of the discussed example CH, similar average velocities were observed 
over all angle space (Figure 4.3[ i). This indicates that there is no preferred 



direction in the boundary evolution of the CH, and hence the boundary evolution 
is isotropic. The value of the average signed velocities, (v), is ~0 km s~^ in 
all angle bins, which indicates no consistent expansion or contraction in the CH 
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360 



Figure 4.3: (a) The CH contour observed at 05:48 UT (solid line) and 17:00 UT 
(dashed line) on 2008 September 8. The frame of reference is the earlier contour 
centroid and angle is measured anti-clockwise from the right direction of the earlier 
contour centroid. (b) The radial location of the boundary plotted as a function of 
position angle (solid: earlier; dashed: later), (c) Boundary displacement velocities 
as a function of position angle, (d) Boundary displacement velocity distributions 
as a function of position angle for all observations of the CH. The dotted lines in 
(a), (b), and (c) mark the boundary locations at 0°, 90°, 180° and 270°. The error 
bars in (d) correspond to ±1 a for each 20° bin. 
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during the time of observation. Furthermore, the majority of ±1 o" values for all 
angle bins are below 2 km s^^ (as is the case for all eight CHs studied here), which 
is in good agreement with the average velocities of random granular motions. 

In order to extract characteristic velocities, all CH observations were com- 
bined over all angles to obtain the unsigned boundary velocity distribution (Fig- 
4.5^ ). From this distribution typical unsigned velocities were found to be 



ure 



0-2 km s~^ 

4.4.2 Magnetic reconnection rate 

Magnetic reconnection between CH open fields and small-scale QS loops (with 



heights of < 3 Mm; Wiegelmann & Solanki 2004) is expected to occur at the 



transition region level. However, the CH boundary displacements are obtained 
from 195 A images that show transition region and the inner coronal heights. 
Here we can make the reasonable assumption that as a small-scale loop reconnects 
with the open magnetic field in the transition region it will consequently cause 
a displacement between open and closed field lines in the inner corona. In fact, 
the displacement at the transition region should be roughly proportional to the 
displacement detected in the EUV images if the surrounding QS has a relatively 
stable magnetic configuration in the corona and the mentioned reconnection does 
not fundamentally change that. We use velocities obtained from FIT 195 A 
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images as we have no information on the displacement velocity at the reconnection 



site. The Alfven speed, fA, was calculated using = Bj^Amm^n^, where B 
is the magnetic field strength and is the electron number density. Typical 
values for the magnetic field (5 = 2 G) and density (rig = 10*^ cm~^) in the lower 



corona were used from Schrijver & van Ballegooijen (2005 ) to calculate the Alfven 



speed, which was found to be ~600 km s ^ . The magnetic reconnection rate was 
determined from 

Mo = v,^lvK, (4.14) 



where fin is the infiow speed (see Figure 4.4), the relative speed of the opposite 
polarity magnetic field lines being brought together prior to reconnection. The 
distribution of displacement velocities was used for t>in and substituted together 



with the Alfven speed into Equation 4.14 to obtain the reconnection rate dis- 



tribution. This distribution is shown in Figure 4.5 where typical reconnection 
rates are in the range 0-3 x 10^^. 

4.4.3 Magnetic diffusion coefficient 

We calculate the magnetic diffusion coefficient for each 5° segment of each CH 



observation using Equation |4.12| by substituting the boundary displacement dis- 
tances in the segment [br) and the time elapsed (5t) between consecutive ob- 
servations. The resulting distribution of diffusion coefficients is shown in Fig- 
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'out 




Figure 4.4: Magnetic reconnection between a closed and an open field line. The 
solid lines show the field lines pre-reconnection, and the dashed lines demonstrate 
the field lines after the reconnection. The arrows indicate the direction of the inflow 
with velocity Vm, and the outflow with velocity Vout- 



ure 



4^, with values typically in the range 0-3 x 10^^ cm^ s ^. Using typical loop 



footpoint distances and characteristic times for reconnection, Fisk Schwadron 



(2001) have determined the diffusion coefficients inside and outside of a CH to 



be 3.5 X 10^'^ cm^ s ^ and 1.6 x 10^^ cm^ s ^, respectively. 



4.4.4 Magnetic field in coronal holes and the quiet Sun 



Establishing the CH boundaries also allowed us to determine the average magnetic 
field strength (<-Blos>) of the CH and the surrounding QS within the "box". For 
the studied CHs, the signed <i?Los> "was found to be 2.3±0.9 G and -2.3±0.5 G, 



while -0.1±0.7 G in the QS surrounding the CH. Figure 4.6 shows the mean 
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Figure 4.5: Frequency distributions of boundary displacement velocities (a), mag- 
netic reconnection rates (b), and diffusion coefficients (c). 



-Blos distributions from all CH observations and the corresponding QS regions. 



The thick/blue and thin/red lines correspond to the predominantly negative and 
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Figure 4.6: The <-Blos> histogram of the negative and positive polarity CHs 
(thick/blue and thin/red line, respectively) and the QS (dashed line). 



positive CHs respectively, and the dashed black line corresponds to the QS. The 
figure shows a clear separation between the negative polarity holes and the QS, 
while there is some mixing between the positive polarity holes and the QS. The 



top plot in Figure [47| shows the unsigned <-Blos> as a function of the CH area. 
The plot shows that the unsigned <-Blos> in larger CHs varies less. However, 
no correlation has been found between the area and the unsigned <i?Los> of the 



CHs (correlation coefficient ~0.23). The bottom plot in Figure 4.7 shows the 
location of the CHs, their polarity, and the range of the corresponding <i?Los> 
values. 
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Figure 4.7: Top: The CH mean -Blos as a function of CH area. Bottom: The 
latitudinal locations of individual CHs and the CH mean -Blos measures in the 
CHs. Different colours/symbols represent different CHs. 



4.5 Discussion and Conclusions 



The relocation of open magnetic field lines through interchange reconnection is 
the most widely accepted method to explain some of the distinctive characteris- 
tics of CHs - namely their rigid rotation, the lack of loops with heights >10 Mm 
inside CHs, and the CH boundary relocations through interaction with QS loops. 
In this study CH boundary locations were determined using an intensity thresh- 
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olding technique described in detail in Chapter [3j This method determines the 
CH boundaries based on observations rather than allowing the observer to choose 
a subjective threshold. A number of CHs were tracked for several days to study 
the changes in the CH boundaries. The CH boundary tracking was carried out 
by temporally coupling the boundary contours, which allowed relative distances 
of boundary movements and hence velocities to be established. In our study the 
displacement velocities were found to be isotropic with unsigned values typically 
0-2 km s~^, which imply reconnection rates of 0-3 x 10^^. Studying the dis- 
tribution of the velocity, magnetic reconnection rate and the diffusion coefficient 



(Figure 4.5 ) the decay was found to follow a power law, which suggests that these 
properties are scale-free. 

These boundary velocities are in good agreement with the isotropic nature and 



~l-2 km s~^ average velocity of random granular motions (Markov, 1991 ). Hence, 
we suggest that the majority of the reconnection at CH boundaries is likely to 
be initiated through random granular motions. Velocities in excess of ~2 km s^^ 
were occasionally observed during our analysis. Such surges in the boundary 
displacement could be explained by consecutive interchange reconnections in the 
case of favourable loop alignments, whereby a number of loops have opposite 
magnetic polarity to the CH flux in their nearest footpoint to the CH. In close 
proximity the CH open field lines could continuously reconnect with the loops 
and advance into the QS, before being halted by loops that face the CH with the 
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same polarity in their nearest footpoint to the CH. 

Interchange reconnection has been previously discussed by |Fisk &: Schwadron 



(2001 ) as a mechanism that allows the transfer of magnetic flux from one point to 



another (see Equation 4.13). We determine typical values of diffusion coefficient 



to be 0-3 X 10 cm s using Equation 4.12 with every individual CH boundary 



displacement measurement (i.e., from every 5 degree segment of each CH time 
pairing). Although the values derived in this study are comparable to those 



determined by Fisk & Schwadron (3.5 x 10^^ cm^ s^^ inside a CH), it is important 
to note that the results are achieved by different methods. Those presented 
here are based on CH boundary tracking using observations and reflect upon 



the diffusion processes happening at the CH boundary, while Fisk &: Schwadron 



provided estimates of the magnetic diffusion in CHs based on typical loop heights, 
foot-point distances and sizes of supergranular cells. 

The lack of systematic relative motion between CHs and the QS on timescales 
of days indicates that the small CHs studied rotate at approximately the local 
photospheric QS rotation rate. In the interchange reconnection model, QS loops 
continuously reconnect with the CH open field and effectively "pass" through the 
CH (an argument also used to explain the rigid rotation of CHs). This suggests 
that over several solar rotations a relative velocity should be observed between 
the CH and the QS due to the differential rotation of the QS. If we assumed a 
CH to rotate uniformly at a velocity equivalent to the QS rotation velocity at the 
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latitudinal location of the CH centroid, we would expect the CH to shift eastward 
at lower latitudes and westward at higher latitudes so as to compensate for the 
differential rotation (and the shearing effect) of the surrounding QS. Although 
no evidence for this has been observed in this study, we suggest it is merely due 
to the fact that the longest observation period was six days, during which only 
a few degrees shift would occur over a CH latitudinal width of < 35° (bounding 
box height). 

Studying the small-scale changes in CH boundaries can greatly expand our 
understanding of not only CHs but the QS as well. The interaction between a 
CH and the QS can shed light on the processes controlling the magnetic field evo- 
lution in the QS and how open field is accumulated (CH birth) and diffused (CH 
decay) . Although our current understanding of exactly how CHs are created and 
destroyed is still limited, it is likely that the random granular motions in the QS 
allow temporary accumulations of open field. Small CH candidates will then grow 
or decay based on the local loop-emergence rate and the specific orientation of 
QS loops near the CH boundary. The decay of a CH can be influenced by nearby 
QS loops but also by processes occurring inside the CH. As well as magnetic 
bipoles emerging within a CH, small loops admitted into the CH during inter- 
change reconnection could accumulate over a CH lifetime and ultimately lead to 
the fragmentation and diffusion of the CH. Bipoles and open flux regions would 
have to be tracked in time to provide observational evidence for the processes 
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governing the birth and the decay of CHs. This could also provide information 
on how small CHs turn into large, stable and long-lived CHs, and what processes 
initiate the decay of a long-lived CHs. 
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Chapter 5 

High-speed solar wind forecasting 



Recurring high-speed solar wind streams observed at 1 AU are known to be as- 
sociated with CHs. However, the relationship between CH properties and in-situ 
HSSW properties needs to be further investigated to aid space weather forecast- 
ing. In this Chapter, we describe the Minor Storm (MIST) algorithm developed 
by the author in order to link CHs to HSSW streams detected at Earth. Over 600 
CH-HSSW pairs were analysed over nine years which allowed us to determine 
the expansion factor of open magnetic flux tubes rooted in CHs. The expansion 
factor was found to be 2-48, which is similar to values previously found using 
theoretical models. The study also revealed the HSSW stream velocities to have 
a strong correlation with the proton temperatures and Kp indices. Furthermore, 
we confirmed the relationship between CH-related HSSW streams and enhanced 
high-energy electron fiuxes measured in the Earth 's radiation belt during the de- 
scending phase of the solar cycle. 
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5.1 CHs and high-speed solar wind 

Unlike flares and CMEs, the prediction of CH related geomagnetic disturbances 
can be forecasted with high certainty. This is due to three main reasons: 

• CHs are observable for long periods of time, and thus a forecast can be 
made several days in advance before the CH related HSSW reaches Earth. 

• The velocity distribution of the CH-related HSSW streams is known, and 
can be used to determine the arrival-time range. 

• CHs are sources of continuous HSSW flows and thus no eruption time needs 
to be determined as in the case of flares and CMEs. 

These characteristics are used in the present research to forecast the arrival time 
of HSSW streams at Earth. The start and end-time range of the HSSW periods 
are determined based on CH location, however, the end of the HSSW period 
at Earth is less straightforward to estimate. Since flux tubes rooted in CHs 



expand super-radially (Section 5.1.4), the HSSW period end-time depends on the 



expansion factor ( Guhathakurta & Fisher, 1998). In Section 5.1.5 we discuss a 



semi-automated method to determine the expansion factor, which is then used in 
the MIST algorithm to determine the time-range for the plasma arrival and end 
at Earth. The method is run over a long period (2000-2009) to obtain CH and 
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HSSW properties and geomagnetic disturbance indices, which are compared to 



establish any relationship (Section 5.1.6) 



5.1.1 Previous forecasting methods 



5.1.1.1 The Robbins method 



The empirical forecasting model of Robbins et al. (2006) relies on the CH detec 



tion method developed by Harvey & Recely (described in Section 3.2.1 [Harvey &: 



Recely 2002). KPVT He 1 10830 A spectrograms and photospheric magnetograms 
were used to determine the location, size and magnetic properties of CHs, while 
the solar wind data was obtained from the OMNIWeb website, provided by the 
US National Space Science Data Center. 

The heliographic CH maps were first divided into 23 longitudinal sectors in a 
±77° longitudinal and ±60° latitudinal range. In each sector the percentage of 
CH coverage was determined, and cross-correlated with the solar wind speed time 
series. First, a model solar wind time series was created by determining the CH 
area percentage in the 14° sectors. These CH area percentages in each sector were 
then rescaled to agree with the observed solar wind speed, and the resulting linear 
scaling coefficients were used to forecast HSSW velocities. The model was used 
on the period 28 May 1992 - 25 Sep 2003, where for the best one-month periods 
the model was found to be accurate to within 10% of the observed solar wind 



143 



5. HIGH-SPEED SOLAR WIND FORECASTING 



measurements. The linear correlation coefficient was determined to be ~0.38 for 
the 11 years studied. 

The KPVT CH maps were also used to predict the IMF polarity at Earth. 
The pixel values in CHs were set to +1 or -1 depending on a positive or nega- 
tive magnetic field strength, respectively, while QS pixels were set to zero. The 
magnetic field was averaged in the sector to determine its overall polarity, and 
then used as a baseline for the polarity forecasting. Using this method the IMF 



polarity was correctly forecast 58% of the time. In their paper, Robbins et al. 
stated that their model had inherent inaccuracies which were likely to be caused 
by the oversimplification of the magnetic field structure associated with CHs and 
the corresponding solar wind streams. They also assumed that the solar wind 
velocity is related linearly to the CH area, whereas it has been shown that at a 



certain scale size the wind velocity is independent of CH area (Poletto & Suess 



2004). 



Despite its partial success, it is important to note that the Robbins et al. 
method provided a fast and simple approach to HSSW forecasting by solely using 
CH maps and magnetograms. 

5.1.1.2 The Vrsnak method 



Vrsnak et al. carried out two studies relating CHs and HSSW streams. The first 



study used CH properties to forecast solar wind parameters, while the second 
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study forecasted the scale of geomagnetic disturbances based on the previous 
work. 

First, the CHs were identified in soft X-ray images from the Soft X-ray Imager 



aboard the GOES-12 spacecraft (SXl; Hill et al., 2005). The CH boundaries were 



determined by setting a fixed intensity threshold of 0.15 DN s ^. For each day 
the fractional CH area was determined in three 20° wide longitudinal sectors. 
The studied longitudinal ranges were [-40°, -20°], [-10°, 10°], and [20°, 40°] from 
the central meridian. The latitudinal range of the analysis was ±45°. 

The time range of the study was 25 January to 5 May 2005, where eleven 
distinct peaks in the solar wind properties were identified and linked to CHs. The 
HSSW properties measured were the velocity (t>), magnetic field strength (B), 
proton number density {n), and the proton temperature (T), using data from the 
ACE/SWEPAM and MAG instruments. The time when a CH had the largest 
area in the sector was compared to the onset time of the solar wind property 
enhancements. The time delays were quantified by applying a cross-correlation 
analysis. This was followed by the time-lagged correlation of the maximum CH 



areas with the HSSW properties as shown in Figure pA] The obtained coefficients 
were then used in a fully automated manner to forecast HSSW properties several 
days in advance. 

The most accurate prediction was obtained for the solar wind velocity, for 
which the average relative difference between the calculated and the observed 
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Figure 5.1: The time-lagged correlations of the CH areas in the central- meridian 
sector and the peak values of a) solar wind velocity; b) proton temperature; c) 
density and d) magnetic field, for 25 January - 5 May 2005. The least-squares fits 
and the corresponding correlation coefficients are shown in the insets. The time 



lags used are shown in the plot headers ( Vrsnak et al. 2007a ) . 



peak values was 10%. The forecast reliability for T, B, and n was found to be 
20, 30, and 40%, respectively. It was also determined that the peak in HSSW 
properties appear in the sequence: n, B, T, and v with time delays of 1, 2, 3, 
and 4 days, respectively, relative to the CH transit through the central meridian 
(CM). 

Note that in this method no magnetograms were used, which leaves room for 
error in determining whether all low intensity regions included in the study were 
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in fact CHs. While a CH linked to HSSW can be accepted as a proof of the CH 
status, it is the western CH boundary position that should be linked to the start 
of the HSSW period rather than the time of the maximum area attained by the 
CH in a CM sector. In addition, several CHs being present on the disk could 
introduce errors in linking CHs to HSSW streams when using only a CM sector 
based analysis. The arbitrarily chosen CH boundary intensity may also mean 
that obtained CH properties could be falsely dependent on the solar cycle, if the 
method was used on long-term analyses. 

The geomagnetic index (Dst) forecasting method was based on the empirical 



correlation between the CH area/location, and the Dst index (Vrsnak et al. 



2007b). Studying the same period as mentioned above, it was found that the 
highest geomagnetic activity occurs 3.2±1.6 days after a low-latitude CH crosses 
the central meridian, and the amplitude of the Dst dip (i.e. enhanced geomagnetic 
activity) is correlated with the CH area and depends on the dominant magnetic 
polarity of the CH. (Note that the area used here is a fractional CH area in the 
sector located at the CM.) The Dst dip amplitude forecasts were found to have 
a rather low accuracy of 30%, and the time of the Dst minimum could only be 
forecasted with an error of 2 days. 
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5.1.2 Parker spiral and the radial plasma propagation 



In Section 1.3.2 we discussed the heliospheric structure of solar wind streams. 
Here, we readdress the hehospheric propagation of the solar wind plasma in order 
to determine its arrival time at Earth. 

Let us consider a static coordinate system with the Earth at its origin. In this 
system the Sun rotates around its axis at the synodic rotation rate (Q), relative to 
Earth. The Sun's rotation speed relative to the stars is called the sidereal angular 
velocity, from which the Earth's orbital rotation rate (i.e. ~360°/365.25 days ~ 
l°/day) is subtracted to get the synodic angular velocity. This is essentially the 
rotation rate at which the Sun appears to rotate from Earth. The synodic angular 



velocity of the Sun has been determined to be ~13.3°/day (Allen, 1973; Beck 



2000| |Wang fc Sheeleyj |1993[ ). 

The equation to determine the arrival time of an ejected plasma blob from 
the Sun can be arrived at in two different ways. The simplest way is to consider 
an observer on Earth facing the Sun in the stationary coordinate system. In this 
system all material leaving the Sun travels radially. Thus, if a region on the Sun 
located at the Sun-Earth line ejects material (i.e. solar wind streaming from a 
CH), the plasma will travel radially a distance of 1 AU to reach Earth. The 
velocity of the plasma can be treated as constant since the acceleration phase 
occurs under 20 which is relatively small compared to the 214 Rq distance 
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between the Sun and Earth (Kojima et a/., 2004). Hence, the time it takes for 



the plasma blob to arrive at the target (Earth) from the time of the observation 
at the CM can be determined from the simple equation 



dti 



dr 



(5.1) 



where dr is the distance the plasma travels (1 AU), and Vsw is the solar wind 
velocity. Note, that if the source region is not at CM, the time it takes for the 
source to rotate to CM needs to be calculated using 



dto 



da 

T7' 



(5.2) 



where da is the angular distance from the source to the CM. Taking this into 
consideration, the total time from observation to arrival can be written as 



, , , dr da 

dt = dti + dt2 = h -TT- 



(5.3) 



This equation can also be derived using Equation 1.24 



r — ro 



{a - tto), 



(5.4) 



where Ur = Vsw, r — ro = dr and a — ao = da = Qdt2. Substituting these notations 
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Figure 5.2: A poleward view of the radial propagation of plasma elements (61, 62- 
63) along the Sun-Earth line. The elements are ejected from different longitudinal 
positions as the Sun rotates (to, and ^3 represent different times of observation). 
O, is the synodic solar angular velocity. 



leads to Equation 5.2 



Figures 5^, 5^ and 5^ demonstrate the connection between the radially 
propagating solar wind elements and the spiral structure of the streaming solar 



wind. In Figure 5.2 plasma elements (61, 62 and 63) are continuously ejected from 
a solar point-source along the Sun-Earth line, as seen from a poleward view. The 
source rotates with the Sun at the synodic rotation rate, hence, the source of the 
plasma elements propagating towards Earth is different at different observation 



times (to, ti and ts). In Figure 5.3 the location of a single plasma element (red 
blob) is followed at consecutive observation times relative to the source point. 
The yellow blobs represent the previous locations of the plasma element relative 



to the source point. With the rotation of the Sun, the source rotates away from 
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Figure 5.3: The propagation of a single plasma element in time along the Sun- 
Earth line. The red blob shows the plasma position at a given observation time. 
The yellow blobs represent the previous locations of the plasma element relative 
to the source location on the Sun. The coloured radial lines show the radial lines 
relative to the source location at given observation times. At ^3 a spiral structure 
appears to have formed relative to the source location on the Sun. 

the Sun-Earth hne. Hence, if an observer was located in the source point, the 
location of plasma element at different times would reconstruct a spiral path. 



Figure 5^ shows the plasma elements continuously being ejected from an 
extended source region on the Sun. The blobs streaming along the Sun-Earth 
line originate from different longitudinal locations on the Sun, and the colours 
denote different ejection times. Note the spiral geometry taking shape relative 
to the source region due to the rotation of the stream source. This figure also 
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Figure 5.4: Multiple plasma ejections from non-point source on the Sun. Different 
colour blobs correspond to different times of ejection. As the source region rotates 
the relative positions of the ejected plasma elements at consecutive observation 
times recreate the well-known heliospheric structure of the solar wind. 



demonstrates the relationship between the latitudinal width (i.e. dcp) of the source 
region and the duration of the plasma stream observed at the target (Earth). 

Although understanding the heliospheric propagation of the solar wind plasma 
is important, we only need to consider the radial propagation of the plasma to 



provide HSSW forecasts based on CH locations (i.e. Equation 5.3). 
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5.1.3 CH locations and the HSSW 

As discussed in the previous section, to determine the HSSW arrival time at 
Earth the location of the source region has to be determined. In this study it is 
assumed that the magnetic field line of HSSW stream to first connect to Earth 
is rooted in the west-most point of a CH. Here, we demonstrate the relationship 
between the CH location on the Sun and the HSSW periods detected at 1 AU. 



Figure 5.5 shows the CH positions at the start, maximum and end of the 
high-speed solar wind period on 14-24 June 2008. The images show that the 
west-most CH boundary was at approximately W50° when the HSSW stream 
arrived at the ACE spacecraft. The CH has passed over the west limb by the 
time the solar wind speed dropped back to normal levels (300-400 km s^^). 

The relationship between the CH and its associated solar wind stream mea- 
sured at three distinct points in the Heliosphere was further explored using imag- 
ing and in-situ measurements from STEREO B, ACE and STEREO A. The top 



row of Figure 5.6 shows the CH complex on three days: 12, 14 and 16 June 2008. 



The corresponding solar wind velocity measurements, given in the bottom 



three panels of Figure |5.6[ show a high-speed solar wind stream sweeping around 
the inner Heliosphere, from STEREO B to ACE to STEREO A. The HSSW 
stream was initially detected by STEREO B at 12 Jun 2008 19:10 UT, when the 
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Figure 5.5: Top: EIT Lambert maps showing the CH position on the day of the 
start, the maximum and the end of the HSSW period. Bottom: The corresponding 
in-situ ACE solar wind speed measurements for the studied period. 

west-most boundary of the CH was at ~W53 from the perspective of STEREO 
B. Subsequently, the solar wind speed began to rise in ACE in-situ measurements 
two days later on 14 Jun 2008 11:00 UT. At this point in time, the CH's west- 
most boundary was at ~W50 as seen by SOHO. The high-speed stream finally 
traversed the point in the Heliosphere monitored by STEREO A on 16 Jun 2008 
6:30 UT, at which time the western boundary of the CH was at ~W43. The com- 
plimentary imaging and in-situ capabilities of SOHO, ACE and the two STEREO 
satellites therefore offer a unique opportunity to study the relationship between 
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Figure 5.6: Top: EIT Lambert maps showing the CH position on three different 
days when the high-speed solar wind reached STEREO B, ACE and STEREO A. 
Bottom: The solar wind profiles obtained from the three satellites, with the line 
marking the start of the high-speed solar wind. 

CHs rotating across the solar disk and high-speed solar wind streams sweeping 
through the Heliosphere. 
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5.1.4 Flux tube expansion 

It is known that magnetic field lines originating from CHs reach interplanetary 
distances. The field lines extend from within the confined boundaries of CHs and 
form a flux tube (or solar wind stream). It has been shown that these flux tubes 
expand super-radially and the expansion factor can be derived from the mass 
conservation equation: 

— = Ap., (5.5) 

where dM/dt is the mass loss rate, A is the cross-sectional area of the flux tube, 
p is the density and v is the solar wind velocity. Given that the plasma flow is 
constant, the mass-loss rate will be constant, and hence 

Apv = const. (5.6) 

This is the one-fluid description of a proton-electron solar wind. To give the 
cross-sectional area of the flux tube at any point along the tube (s), the equation 



A = Aof{s)[^y (5.7) 



is used, where sq and s are the radial distances where the area is measured, and 



/ is the non-radial flux tube expansion factor ( Guhathakurta et ai, 1999b). 
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Equation |5.7| is modified to contain the observed parameters, namely the 
latitudinal width of the solar wind source and the HSSW duration measured at 
Earth. This is done by assuming the flux tube to have a circular cross-section 
and denoting the 

Ai = nrl (5.8) 
for the radially expanded flux tube, and 

A2 = Tirl (5.9) 

for the superradially expanded flux tube. Here, ri and r2 are the flux tube radii for 
the radially and superradially expanded flux tubes, while Ai and A2 are the flux 
tube cross-sectional areas for the radially and superradially expanded flux tubes. 
If we assume a radial expansion, with the flux tube having a A0i latitudinal 
width at Earth, the width can be approximated as sA0i = 2ri. Substituting ri 



into Equation |5.8| hence yields 

7rs2(A0 



A, = HIA^. (5.10) 



This is substituted into Equation 5.7, for the radial expansion (/ = 1) to get 



^£W.^. (5.n) 

4 so^ 
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Applying the same steps, by substituting Equation |5.9| and sA(f)2 = 2r2 into 



Equation 5.7, for superradial expansion (/ 7^ 1) we get 



(5.12) 



These latter two equations can be rearranged to yield 



A01 



4Ao 



A02 = a//(s 



4Ao 



(5.13) 



which can then be related as 



(5.14) 



If Ar2 is the HSSW duration measured at Earth, we can assume A02 = ujAt2, 



which we substitute it into Equation 5.14 to express the expansion factor 



/ 



U)At2 

A0r 



(5.15) 



As a demonstration. Figure 5.7 shows the Sun (yellow disk) with a solar wind 
source of A0i latitudinal width, and the expanded flux tube with a A02 width 
at Earth. The solar wind is shown to expand radially and is measured to last for 
Ati at Earth (green disk), while the super radial expansion is shown to give rise 
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Figure 5.7: Schematic showing the super-radial expansion of a CH flux tube. 
The CH has a latitudinal width A(/>i, the radial expansion is shown with the inner 
black lines extending from the CH, the the outer red lines indicate the super-radial 
expansion. 

to a longer solar wind period of Ar2 at Earth. 

5.1.5 Expansion factor: semi-automated approach 

In order to determine the HSSW stream expansion factor, the properties of 47 
CHs were studied in the period 2005-2008. The CHs were identified and grouped 
using the CHARM algorithm, and the west-most boundary locations were used 
to determine estimated HSSW arrival times at Earth assuming an average solar 
wind velocity of 600 km s~^. 
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Figure 5.8: Scatterplot showing the relationship between latitudinal width of CHs 
and the HSSW duration. 



Figure 5.8 shows the correlation between the CH latitudinal widths and the 



corresponding HSSW durations of the 47 CHs analysed in this study. The green 
line was fitted to the data point with the largest gradient, while the red line was 
fitted to data point with the smallest gradient. 

An estimated end-time was also determined using the east-most CH boundary 
location, assuming a simple radial expansion. The estimated arrival and end-time 
was used to avoid confusion between multiple HSSW periods in the time range. 
With the help of the estimated arrival and end-times, the HSSW period start and 
end was manually determined using solar wind velocity data obtained from the 



ACE/SWEPAM instrument. Figure 5.9 shows an example of a CH automatically 
detected on 13 March 2007, and the corresponding HSSW period identified in the 
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Figure 5.9: The automatically detected CH group (SOHO/EIT 195 A Lam- 
bert map) and the semi-automatically identified corresponding HSSW stream 
(ACE/SWEPAM). 



ACE/SWEPAM solar wind data. The dashed red lines show the estimated start 
and end of the HSSW period (left and right respectively), based on a 600 km s~^ 
solar wind velocity and a radial expansion. The green lines show the detected 
HSSW start and end time. 
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Figure 5.10: Histogram showing the frequency distribution of the expansion fac- 
tors obtained from 47 CHs. 



Using Equation 5.7 and the line fit gradients (Ar/A0) shown in Figure 5.8, the 



expansion factor was determined to be in the range 5-48. Figure 5.10 shows the 
frequency distribution of the expansion factors, where the mean expansion factor 
was found to be 15. However, the most frequent expansion factors were found to 
be <10, while large expansion factors (>10) were found to be rare in this study. 
The expansion factor range determined in our study was found to be similar 
to previous theoretical and observational results that find the expansion factor 



to be up to 20 at lAU DeForest et al. (2001); Guhathakurta & Fisher (1998); 



Guhathakurta et al. 


( 


1999b 


); 


Wang 


( 


1994 


) . Figure 


5.11 
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Figure 5.11: CH magnetic flux tube expansion factors as a function of distance 



for ten co-latitudes from 0° to 30° ( Guhathakurta et al, 1999b). 



factors inferred by Guhathakurta et al. (1999b) from Ulysses observations at 1 
AU. 



5.1.6 Expansion factor: automated approach 

Using the expansion factor determined from the semi-automated approach, a 
similar analysis was carried out in an automated manner to obtain the expan- 
sion factors using a larger number of CHs and HSSW observations. Using the 
west-most boundary position of a CH, and a velocity range 400-900 km s~^, an 
approximate arrival time range was determined, where the rising phase in the 
HSSW velocity was automatically identified. Similarly, the estimated end-time 
range was determined from the east-most CH boundary location and previously 
obtained expansion factors. In this end-time range the HSSW end is determined 
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where the velocity reaches minimum. Figure |5.12| shows an example CH group 
which was identified on the 13 January 2007 in an SOHO/EIT 195 A image, with 
the corresponding HSSW stream identified in ACE solar wind velocity data. Note 
the very similar onset times of the v, T, p and Kp enhancement. 

In the present analysis the following criteria were used: CH groups have to 
be of substantial size (> 1.5x10^ Mm), the HSSW rise phase has to occur within 
~2-4 days of the CH boundary appearing at the CM (calculated from suggested 
SW velocities of 400-900 km s^^), the HSSW velocity enhancement has to be 
at least 150 km s~^, and the maximum velocity has to be above 400 km s~^. 



An example is shown in Figure 5.12, where a large CH group is shown with the 
corresponding HSSW stream detected in ACE solar wind data. In the solar wind 
velocity plots the solid red line shows the approximate HSSW arrival time if the 
solar wind speed is of 600 km s~^. The dashed lines correspond to HSSW arrival 
times for velocities of 900 km s^^ and 400 km s~^ (left and right, respectively). 
The green solid line represents the observed ("true") start of the HSSW period. 
This time was automatically detected based on the previously mentioned criteria. 
The estimated end-time range are indicated with dashed blue lines and the cyan 
solid line shows the automatically detected end-time. The label "CH Group: 1" 
indicates that the detected HSSW period is linked to CH group 1, shown in the 
FIT Lambert map above. 

In the period 2000-2009, 616 CHs were automatically linked to HSSW streams. 
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Figure 5.12: The automatically detected CH group (SOHO/EIT 195 A Lambert 
map) and the corresponding HSSW stream (ACE/SWEPAM bulk velocity). 
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CH latitudinal width [°] 

Figure 5.13: The automatically obtained CH latitudinal widths as a function of 
the corresponding HSSW durations. The green and red line fits are obtained from 
the automated data, the dashed lines corresponding the semi-automated method 
are shown for comparison. 



From the automatically obtained CH and HSSW properties, the new expansion 



factor range was determined to be 2-48. Figure 5.13 shows the HSSW durations 
as a function of the CH latitudinal widths of the 616 CHs studied. The line 
fits from the semi-automated method are shown for comparison in dashed red 
and green lines, while the line fits for the automated method are shown in solid 



lines. Figure 5.14 shows the frequency distribution of the expansion factors, where 
the mean expansion factor was found to be 20. Similar to the semi-automated 
approach, the most frequent expansion factors were found to be <10. It is also 
clear from the histogram that small expansion factors are much more frequent 
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Figure 5.14: Histogram showing the frequency distribution of the expansion fac- 
tors automatically obtained from 616 CHs. 



than large expansion factors. 



5.2 Results 



5.2.1 CH area and HSSW properties 



Previous work by (Vrsnak et ai, 2007a) has shown HSSW properties to be in- 



fluenced by the size of the corresponding CH (see Figure 5.1 in Section 5.1.1.2) 



In the study of [Vrsnak et al.\ CH areas were determined in a CM sector of 20° 
in the time period 25 January to 5 May 2005 and were compared to peak values 
of HSSW properties. In our study CH areas were determined for CH groups 
without sectorial limits in the time period 10 August 2000 to 24 June 2009. We 
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found a weak correlation of ~0.15 between the CH group area and the maximum 
HSSW velocity. No correlation was found between the CH area and HSSW rimax, 



\B 



z,max I ; 



Kpmax, Tmax- These results are in disagreement with those acquired by 



(Vrsnak et ai, 2007a). The scatter plots and the frequency histograms of the 



HSSW properties are shown in Figures 5.15, 5.16 and 5.17 



The correlation between the CH area and the mean HSSW properties was 
also studied. A weak correlation of ~0.2 was found between the area and <T>, 
<Kp> and <v>. No correlation was found between the area and <n>, <|i?^|> 



(Figures [5l8| [57l9| and [Oo| . 

When inspecting the relationship between the CH area and HSSW duration, a 
strong correlation was found (correlation coefficient of ~0.58); this is in agreement 



with the super-radial flux tube expansion model (Figure 5.20) 



When relating the maximum HSSW velocity to the maximum temperature 
a strong correlation of ~0.61 was found, and and even stronger correlation was 
shown between the mean HSSW velocity and mean temperature (correlation coef- 
ficient of ~0.89). This correlation may indicate that the mechanism accelerating 
the solar wind in CHs is also responsible for the heating of the plasma (e.g. Alfven 
wave heating). 

The mean HSSW velocity and the mean Kp index have a correlation of ~0.49. 
The rest of the HSSW properties show little or no correlation. The described 



scatter plots are shown in Figure 5.21 The frequency histograms of the discussed 
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properties are shown in the corresponding figures. 
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5.2.2 Geomagnetic effects of HSSW streams 

The CH location, area and magnetic flux are known to follow the solar cycle 



(Wang et a/., 1996). Hence, it is safe to assume that the properties of HSSW 
streams originating from CHs should consequently show some variation over the 
solar cycle. The space weather effects of CHs and the corresponding HSSW 
streams is particularly important in this study. A key aspect of the analysis 
presented in this section is the variation in the Van Allen radiation belts, which 
have been linked to solar storms. Relativistic ions and electrons trapped in the 
Earth's magnetosphere make up the Van Allen belts, which have been shown to 



change under the influence of magnetic and electric fields. Figure 5.22 shows 
the radiation belts, which are made up of two layers - the inner (~1.5 REarth) 
and outer (~4-7 REarth) belt. Geomagnetic storms caused by CMEs and CIRs 
cause enhancements in the radiation belt relativistic electron flux, which can 



cause serious damage to electrical devices on satellites (Miyoshi & Kataoka, 2011 



Pilipenko et al. , 2006). It has been shown that during solar maximum activity the 
anomalies experienced on satellites are mostly caused by energetic solar protons 
originating in CMEs, while in the declining and growing phase of the cycle it is the 
magnetospheric electrons that cause satellite disturbances. The outer radiation 
belt has been found to be influenced mainly by CIRs, while the inner belt by 
CMEs. 
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Figure 5.22: The inner and outer Van Allen radiation belts. 

In this study we show the variation in the HSSW properties over nine years 
from 2000 to 2009, which has shown an interesting trend in the CH area, the 
mean HSSW velocity, mean temperature and mean Kp index. Near the start 
of 2003 all four properties show a remarkable rise (hence the relatively strong 



correlation described in the previous section). This is shown in Figures 5.23 and 

EM 

Similar enhancements have been shown in the yearly averaged > 300 keV elec- 



tron flux measured in the Earth's Van Allen radiation belt (Miyoshi & Kataoka 



2011). Figure 5.25 compares the variation in the sunspot number (blue line) dur- 
ing the solar cycle and the electron flux (red flux) measured in the outer portions 
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of the outer radiation belt by the MEPED instrument aboard the NOAA/POES 
satellites. The radiation belt shells are described with "L", which corresponds 
to a set of magnetic field lines which cross the Earth's magnetic equator at a 
distance equal to the L-value measured in Earth-radii. In all outer radiation belt 
shells (especially L=5 and 6) the high energy electron flux rises during the declin- 
ing phase of the solar cycle. The maximum value is reached between 2003-2004 
which is in perfect agreement with the enhancements shown in the CH area, the 
corresponding HSSW stream velocity, temperature and the Kp index of the re- 



sulting geomagnetic storms (Figure 5.23). Although the rise is not as significant 



it can still be recognised in the HSSW density and \Bz\ (Figure 5.24). The CH 
area, the mean HSSW velocity and the temperature show a significant rise which 
is followed by a similarly steep decline to the end of 2005. The average Kp index 
however stays at ~3 from 2003-2006. 

The peak values of the HSSW streams were also studied over the nine years, 
but no significant rise was identified in the year 2003. However, the maximum 
Kp index shows a significant enhancement from 2003 to 2005, reaching a peak 
value of Kp=7.5 in 2005-2006, which is followed by a steep decline at the start 



of 2006, after which it levels off at Kp ~4 (Figure 5.24) 



Figure 5.26 of Miyoshi & Kataoka (2011) shows the geomagnetic storm Kp 



indices corresponding to different causes from 1996-2008. Panel (a) shows the 
frequency of CME related storms, with a small peak in 2002, which coincides with 
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the solar cycle maximum. Larger CME related storms appear to peak between 
2000-2001, also the time of the maximum solar activity - as shown in panel (b). 
Panel (c) shows the CIR related storm frequency enhancement from 2002-2006, 
which is in agreement with the CH related maximum Kp indices found in our 
study. Panel (d) and (e) show the CH related geomagnetic storm frequency 
with significant enhancements in 2002 and 2003, which were also found to be in 
agreement with our results. 
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Figure 5.25: The yearly average > 300 keV electron flux (red line) from 1979- 
2008 for different layers of the Earth's radiation belt (L=4, 5 and 6). The sunspot 
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numbers are shown with a blue line for comparison (Miyoshi & Kataoka, 2011). 



5.3 Discussion and conclusions 

CIRs created by the HSSW streams originating in CHs are well known to cause 
long-lasting geomagnetic storms due to the high plasma velocity and magnetic 



field fluctuations (Borovsky & Denton, 2006). While the geomagnetic storms 
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Figure 5.26: The frequency of geomagnetic storms corresponding to different 
sources: CMEs (a & b), CIRs (c), and CHs (d & e, for north and south directed 



magnetic fields, respectively). Figure courtesy of (Miyoshi & Kataoka, 2011). 



associated with CHs are normally of a minor scale {Kp=4-5), the number of 
occurrences and duration of these minor storms make them comparable to CMEs 
in their contribution to the overall geomagnetic activity. Minor storms are espe- 
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cially important to consider during the declining phase of the solar cycle when 
the number of their occurrences and the scale of the disturbances rise. 

In order to automatically detect and forecast CH related HSSW streams, 
the author of the present thesis developed the MIST algorithm, using which 
the expansion of HSSW flux tubes was estimated. The CH latitudinal widths 
were used in conjunction with the HSSW durations to calculate the super-radial 
expansion factor which was found to be 2-48. Although this result is similar to 



the values of 2-20 determined previously using theoretical models (Wang, 1994), 
it is yet to be investigated why the observational values give a much higher upper 
limit. However, it is likely that the circular fiux-tube approximation based on the 
maximum latitudinal width of a CH will introduce some variation in the results, 
since CHs are known to often have complicated morphologies. Furthermore, in 
our analysis HSSW streams are linked to CH groups rather than individual CHs, 
and hence the latitudinal width of the stream is related to the width of the CH 
group. 

The relationship between CHs and the HSSW streams and geomagnetic dis- 
turbances was also investigated. For the first time, a large-scale study involving 
nine years of data from 2000-2009 linked CH occurrences to geomagnetic storms 
and radiation belt electron flux enhancements. An enhancement was found in 
the CH area, mean HSSW velocity, temperature and Kp index in the year 2003, 
the same time when the high-energy electron flux measured in the Earth's outer 
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radiation belt rose significantly. It was also shown that while the average Kp 
indices reached a maximum value in 2003 and plateaued until 2006, the most 
geoeffective storms reaching Kp = 7.5 occurred in 2005, and kept occurring for 
approximately one year. It is yet to be investigated what causes the enhanced 
geoeffectivity associated with CHs during the declining phase of the solar cycle, 
and which CHs are likely to cause such strong storms. 

Forecasting HSSW streams related to CHs is also beneficial for the calculation 
of the interplanetary CME (or ICME) arrival times at Earth. It is well known 
that CMEs are slowed down as they pass through the ambient slow solar wind 
(also known as the "snow-plough" effect), while this deceleration reduces when 
a CME passes through the high-speed, low density plasma originating from CHs 



Odstrcil et al. (2005); Tappin (2006). The morphology of the CME front could 
also be modified due to interaction with the CIR. For these reasons, it is very 
important to continue research on CH related HSSW streams. 

The author would like to thank Edward Erwin (NOAA/NGDC) for providing 
the Kp data used in this study. 
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Conclusions and Future Work 



The goal of the research presented in this thesis was to study the mechanisms con- 
trolling coronal hole evolution, and to determine the space weather effects of the 
high-speed solar wind streams originating from coronal holes. In order to carry 
out the research, the author developed an automated detection method to stan- 
dardise coronal hole identification throughout the solar cycle and allow the use of 
different instruments and wavelengths. It consequently allowed the investigation 
of the interaction between coronal holes and the quiet Sun, which reflects on the 
mechanisms controlling the short-term evolution of coronal hole boundaries. The 
study of CHs was extended to analyse the long-term changes in coronal hole and 
corresponding high-speed solar wind stream properties. In addition, a method was 
developed to provide forecasts of coronal hole related geomagnetic storms several 
days in advance of their arrival at Earth. In this Chapter, the primary results of 
the analyses are summarised and discussed. Suggestions for improvements and 
future work are addressed at the end of the Chapter. 



187 



6. CONCLUSIONS AND FUTURE WORK 



6.1 Thesis summary 

The goal of the research presented in this thesis was to investigate the behaviour 
of CHs from their small-scale boundary motions seen in solar coronal images to 
the geomagnetic effects experienced at Earth due to the HSSW emanating from 
them. 

It was shown that the fast and robust detection method ( CHARM) developed 
by the author allows the objective identification of CHs throughout the solar cycle. 
Using an additional package {CHE VOL) we were able to study the boundary dis- 
placements in CHs and determine the diffusion coefficient (fi;=0-3xl0^^ cm^ 
as well as the magnetic reconnection rate (Mo=0-3xl0~^). The results indicated 
that the interchange reconnection between CHs and QS loops is likely to be ac- 
commodated by random photosphcric granular motions. Our results support the 
interchange reconnection model at CH boundaries. 

The MIST algorithm allowed CHs to be linked to HSSW streams in order 
to study their geomagnetic effects. An extensive study was carried out, which 
included over 600 CH and HSSW stream pairs. Using the CH latitudinal widths 
and HSSW durations measured at Earth we were able to determine the open 
flux tube expansion factor (/=2-48), which was found to be similar to those 
derived from theoretical studies. The CH and HSSW properties were compared 
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to determine any relationsliip, however, the CH area was only found to correlate 
with the duration of the HSSW streams. Nevertheless, a strong correlation was 
found between the solar wind velocity and temperature, which could support the 
theory of Alfven wave heating in CHs. 

Furthermore, long-term trends were discovered in the CH related HSSW 
stream properties, which were linked to enhanced electron fluxes measured in 
the outer Van Allen radiation belt. This confirms that the the Earth's radiation 
belt is strongly affected by the HSSW streams originating from CHs, especially 
during the descending phase of the solar cycle when the frequency and scale of 
CH related disturbances increases. 



6.1.1 Principal results 

The aim of the work presented in this thesis was to further our understanding 
of the short-term and long-term evolution of CHs. The small-scale changes were 
investigated to reflect on the mechanisms controlling CH growth and decay, while 
the long-term study showed the relationship between CHs and the corresponding 
HSSW streams and geomagnetic effects. Here, we summarise the results obtained 
from these analyses. 

• The Ci/74i?M algorithm was developed by the author to provide an objective 



approach to CH detection throughout the solar cycle (Krista & Gallagher 
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2009). As opposed to previous methods where arbitrary thresholds were 
chosen, here, a local intensity thresholding method was used to determine 
the CH boundary intensity. Additionally, the flux imbalance was deter- 
mined using magnetograms to establish whether the detected region was a 
CH. The detected CHs were then grouped to aid the forecasting algorithm 
and to help the future tracking of CHs. As part of the algorithm the CH 
properties (location, boundary positions, area, line-of-sight magnetic fleld 
strength and intensity) were stored for analysis. The developed algorithm is 
fast and robust, and allows the use of different instruments and wavelengths. 

The CHARM algorithm was used in conjunction with the CHE VOL pack- 
age to study the short term changes in CH boundaries. The study aimed to 
compare the observational results to those derived from the interchange re- 
connection model. The model describes the "relocation" of open magnetic 
field lines through reconnection between open and closed field lines. This 
mechanism can explain the observed changes occurring in the CH bound- 
aries on timescales of hours. 

The results show that the CH boundary displacements were isotropic and 
had velocities in the range 0-2 km s~^. This is in agreement with the 
velocity and isotropic nature of photospheric granular motions. The author 
suggests that the interchange reconnection is aided by the random granular 
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motions and allows small-scale loops get closer to open field lines, and 
accommodate interchange magnetic reconnection. This mechanism could 
also play an important role in the birth, evolution and decay of CHs. It is 
still unclear what initiates the formation of CHs and what determines the 
location of CH formation, however, it is possible that random concentrations 
of open field lines are aided by granular motions. Similarly, random granular 
motions could also be responsible for the fragmentation and diffusion of 
open field in decaying CHs. 

CH boundary relocation speeds in excess of 2 km s~^ were also observed 
which may indicate that reconnection may be accommodated by mecha- 
nisms other than the granular motions. A bipole emergence near a CH 
boundary may cause continuous reconnection between the open and closed 
field lines at a faster rate, while favourable loop alignments could also lead 
to a faster rate of boundary displacements. However, to confirm these the- 
ories it is necessary to analyse images of higher-temporal cadence than the 
12 minute EIT images. 

Using the interchange reconnection model the diffusion coefficient was cal- 
culated to be 0-3x10^^ cm^ s^^ at CH boundaries. This was found to 



be in agreement with the theoretical results derived by Fisk (2005), who 



determined the diffusion coefficients inside and outside of a CH to be 
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3.5 X 10^^ cm^ s ^ and 1.6 x 10^^ cm^ s ^, respectively. 

Using the CH boundary displacement velocities and the Alfen speed typical 
for the lower corona, the magnetic reconnection rate (or magnetic Mach 
number) at CH boundaries was determined to be in the range 0-3 x 10^^. 
In comparison, a reconnection rate of 0.001-0.01 was found by |Isobe et al. 



(2002) in the case of a long duration event flare in a cusp-shaped arcade. 



As shown in Figure |6.1[ the estimated coronal magnetic field strength is 
similar to that used in our study (~10 G), however, the inflow and Alfen 
speeds are significantly higher in the case of the flare. 

The MIST algorithm was developed to build on the CHARM algorithm. 
It allows CHs to be linked to HSSW streams, and it can also be used to 
provide HSSW forecasts several days in advance. Using this algorithm over 
600 CHs were linked with HSSW periods between 2000-2009. The deter- 
mined CH latitudinal widths and the HSSW period durations were used to 
calculate the expansion factor of the open flux tubes rooted in CHs. This 
was determined by the author to be 2-48, while previous theoretical results 
indicated a range of 2-20. It is not yet clear what causes the difference 
between the results, but it may be due to the approximations used in our 
study (assuming a circular morphology for CHs and the related magnetic 
flux tubes). 
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Figure 6.1: The coronal magnetic field strength, the inflow and Alfen speeds, and 
the magnetic reconnection rates determined for the long duration event flare of 12 



May 1997 (Isobe et al, 2002) 



The properties of the CHs and the corresponding HSSW periods identified 
between 2000-2009 were compared to establish any relationship in order to 
aid future forecasting. However, no correlation was found between the CH 
area and the HSSW properties, except for the HSSW duration. Neverthe- 
less, a very strong correlation was found between the mean HSSW velocity 
and temperature. This could confirm that the mechanism accelerating the 
solar wind in CHs is also responsible for the heating of the plasma (e.g. 



Alfven wave heating). It has also been shown that HSSW streams of higher 
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mean velocities are more likely to cause geomagnetic storms with higher 
Kp indices (correlation coefficient ~0.5). The correlation found between 
the CH area and the duration of the HSSW stream experienced at Earth 
can also aid future forecasting purposes (at present the HSSW duration at 
Earth is determined based on the latitudinal CH width and the expansion 
factor) . 

• The long-term analysis of CHs and HSSW properties allowed the deter- 
mination of solar cycle trends in the properties of the CH related HSSW 
streams. The CH area, mean HSSW velocity, mean temperature and mean 
Kp index were all found to rise in the descending phase of the solar activity 
cycle. The peaks were reached within the start and middle of 2003, and the 
values reduced significantly in the consecutive two years. 

• The maximum Kp index showed a significant rise in 2003, reaching a peak 
value of Kp—7.5 between 2006-2007. This was followed by a steep decline 
in 2006-2007 where is levelled off at Kp ~4. Note that a geomagnetic 
storm of Kp—7.5 is considered "strong", causing power fluctuations and 
long lasting storms (typical of CHs) may result in transformer damage on 
electrical grids. Spacecraft operations are also likely to be affected during 
such storms. 

• The time of elevated HSSW properties were found to coincide with the time 
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of enhanced high-energy electron flux measured in the outer region of the 



Earth's Van Allen radiation belt. Miyoshi & Kataoka (2011) previously 



suggested that CHs may be responsible for this effect, but no systematic 
long-term analysis has been carried out before the study presented in this 
thesis. Based on the author's analysis it is clear that during the descending 
phase of the solar cycle CHs are responsible for the electron surges in the 
radiation belt, which can also cause serious damage to electrical devices on 
satellites. 



6.2 Future work and outstanding questions 
6.2.1 CH detection 

One of the main goals of this thesis is to improve previous CH detection meth- 
ods by determining the CH boundary threshold using local intensity histograms. 
However, it is important to note that similarly to previous methods a global in- 
tensity threshold is used in a single image, while CH boundary intensities could be 
different from each other. Furthermore, different parts of a single CH boundary 
could be of different intensities too. Thus, individual CH boundary thresholds 
could be an improvement to the present detection method. 

For further improvement of the algorithm, it would also be beneficial to test 
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other segmentation methods which could potentially improve or supplement the 
local intensity thresholding, i.e., the watershed operation can determine bound- 
aries between areas of similar intensities based on points where single pixels sepa- 



rate similar areas ( Gonzalez & Woods , 2002 ) . This method proved to be successful 



on single CH study cases, but needs improvement in order to detect CHs through 
the solar cycle. Optimal thresholding could be another approach, where a mul- 
timodal intensity distribution of an image is fitted with Gaussian distributions 



to determine thresholds between areas of different intensity ( Harra et al. , 2000 ) . 
Using this method it is important to a priori determine the number of distri- 
butions to be fitted and the approximate locations of the distributions. These 
properties will have to depend on the time of the solar cycle and the number 
of different intensity phenomena present on the Sun that correspond to different 
intensity distributions. Laplacian filtering could also be a useful approach to de- 
tect CH boundaries. By computing the second derivatives of the image it helps 
to identify whether the change in adjacent pixel values is due to an edge or a con- 
tinuous progression. Applying this approach to the intensity histograms rather 



then an image could also prove useful to determine boundary intensities (O'Shea 



et a/., 2000). Multiscale analysis of solar images could also be used in order to 
detect CH boundaries, by identifying edges at different scales and excluding those 
corresponding to noise. 

CHs are known to be dominated by a single polarity and for this reason 
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determining the CH status of a low intensity region is dependent on the flux 
imbalance. At present the CHARM algorithm uses an arbitrarily chosen value as 
a threshold. Further analysis needs to be carried out to justify this value and to 
possibly correct it. The CH grouping algorithm could also be further developed, 
as at the present the grouping is based solely on the distance between individual 
CHs. It may be advisable to group CHs based on their polarity as well as their 



distance, since regions of similar polarity may be interconnected (Antiochos et al. 



2007) and could be treated as a group from a physical perspective. Another 
addition to the CH detection and grouping algorithm could be the identification 
of single CHs (or groups) in order to study their evolution over several rotations. 
The study of individual CHs could reveal what locations CHs are more likely to 
form at, and what determines the maximum size and the lifetime of a CH. It 
could also prove useful in determining the long term geomagnetic effects of single 
CHs (or CH groups) over several rotations. 

It would also be very beneficial to carry out a long-term CH study using X-ray 
observations and compare the obtained CH boundaries with those determined 
from EUV observations. With the advances of solar imaging, it will be very 
important to study CHs using narrow-band filters that correspond to different 
distinct temperatures, as this could reflect on the three-dimensional morphology 
of CHs in the inner corona. The 3D morphology could also be investigated by 
observing CHs from different spatial angles, i.e. using SOHO as well as STEREO 
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A and B observations. In addition, it would be very useful to carry out a long- 
term study on the temperature, density and and element abundances of individual 
CHs using spectral data and compare the CH properties to corresponding in-situ 
HSSW measurements. 

6.2.2 Short and long-term CH evolution 

A further improvement of the CHE VOL algorithm could be the use of high- 
cadence EUV (or X-ray) imagery. The recently launched SDO mission provides 
high temporal and spatial cadence EUV images that will be very useful in the 
future not only for short-term studies but also for long-term studies of CHs, once 
a large database is built up. 

The short-term CH evolution study raised the question of bipole evolution. 
Since bipoles influence CH evolution through interchange reconnection, it would 
be useful to track the emergence and the displacement of individual bipoles and 
to determine their effect on the CH boundary. It would also be important to 
determine whether bipoles "pass through" the rigid rotating CH, as they move 
with the differentially rotating photosphere, or whether they remain inside the 
CH once admitted through interchange reconnection. This could also reflect on 
the role of bipoles in the decay of CHs. 

As part of the long-term analysis, the author investigated the latitudinal dis- 
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Figure 6.2: The polarity-wise separation of CHs over a solar cycle. The latitudinal 
location of CH centerpoints are shown as a function of time between 2006-2009. 
Red symbols correspond to positive polarity, while blue corresponds to negative 
polarity CHs. 



tribution and polarity of CHs from 1996-2009. Figure 6^ shows a very clear 
polarity- wise separation in CHs from ~ 1996-1999 and ~2005-2009, the ascend- 
ing and descending phase of the solar cycle, respectively. During these phases 
the solar hemispheres appear to be dominated by CHs of the same polarity as 
the polarity of the corresponding pole. This clearly dipolar configuration of the 
global solar magnetic field breaks down during the solar cycle maximum, when 



the global magnetic field becomes multipolar (Bilenko 2002). This can be seen 



in Figure 6.2 between 1999-2005, where the latitudinal location of CHs is mixed 
in both hemispheres. 

It is also interesting that while CHs can form in the vicinity of decaying 
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ARs with diffusing unipolar magnetic fields, ARs are not necessarily the reason 



for CH occurrence (Bilenko 2004). The dynamics of non-AR CHs is likely to be 
determined by the global magnetic field. We suggest further analysis to be carried 
out on the polarity-wise separation of CHs, which may reflect on the mechanisms 
involved in the dipole reversal. 

6.2.3 Space weather 

The MIST algorithm uses the simplified approach of considering the maximum 
latitudinal width of CHs to determine the HSSW duration, without considering 
the latitudinal location of CHs. It is a challenge to determine the shape of the 
open flux tube at Earth and to forecast at what part of the flux tube the Earth 
is likely to pass through. However, it is commonly known that open field lines 
expand towards the ecliptic, and thus the method used is a good approximation. 
However, it is important to note, that the location of the equatorial current sheet 
can significantly modify the location of the flux tube at 1 AU and hence cause 
altered HSSW durations. 

The equation used to calculate the expansion factor is based on the CH area 



(Equation 5.7). This equation was rewritten to include the latitudinal width of 



CHs and the duration of HSSW streams at Earth (Equation 5.15). In doing this 



we approximated the CHs to have a circular shape and their maximum latitudinal 
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width to be related to a flux tube of a similarly circular cross-section. However, 
it is clear from observations that CHs can have elongated shapes and an intricate 
morphology. For this reason it would be useful to further develop this calculation 
to account for the non-circular shape of CHs, which may lead to more reliable 
expansion factors. 

The forecasting algorithm at present uses several criteria (i.e. CH area, HSSW 
velocity difference from start to peak value) that limits the number of CHs and 
HSSW streams to be linked. With further testing these criteria could perhaps be 
reduced or even relaxed to allow a larger number of CHs to be included in future 
studies. 

It would also be very useful to extend the analysis on an even longer time 
period and determine whether the CH related HSSW properties show the same 
enhancements in the descending phase of the solar cycle as shown in 2003. 

The forecasting accuracy will also need to be compared to those given by 
theoretical models. The Community Coordinated Modeling Centei]^ (CCMC) is 
a multi-agency partnership that allows access to solar, heliospheric, ionospheric 
and magnetospheric models. The heliospheric models may also be used for fore- 
casting purposes. Model requests can be submitted to CCMC to map the ambient 
solar corona and inner heliosphere for a Carrington Rotation selected by the user. 

These models provide information on the solar wind velocity, temperature, den- 
^ http://ccmc.gsfc.nasa.gov/index.php 
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Figure 6.3: The ENLIL-WSA model of the solar wind velocity in the heliosphere 
from a poleward view (left), ecliptic cone- view of ±45° from the Sun (middle), and 
the full ecliptic view (right). The white line represents the heliospheric current 
sheet location, and the Earth location is shown with a yellow circle. The date on 
the right is the start time of the model run, the date on the left is the time stamp 
of the diagrams shown. 



sity, and magnetic field strength in the heliosphere, from a poleward view as well 



as a cone- view. 



Figure 6.3 shows the result of the ENLIL-cone model run request. In the left, 
the radial solar wind velocity is shown in the heliosphere from a poleward view. 
In the middle, the cone-model shows the heliosphere in the Sun-Earth direction 
in an ecliptic view of ±45° from the Sun. The diagram in the left shows the full 
ecliptic view of the heliosphere (i.e. a 360° extension of the middle figure). In all 
diagrams the white line represents the location of the heliospheric current sheet. 
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and the location of Earth is shown with a yellow circle. 

ENLIL uses the Wang-Sheeley-Arge (WSA) model which extends a potential 
field source surface (PFSS) magnetic field to 21.5 Rq, past the critical point 
where the plasma velocity exceeds the sound speed. It uses a heliospheric current 
sheet model and a slow and high speed solar wind distribution depending on 
CH locations which are determined from magnetograms. ENLIL results map the 
heliosphere from 21.5 Rq - 1.6 AU, and from -58 to +58 ° heliographic latitude. 

The CCMC heliospheric models all use Carrington rotation magnetograms 
to identify open and closed field regions on the Sun in order to carry out the 
magnetic field extrapolation. This however means that it is an approximation 
of the heliospheric magnetic field, since it uses data which is time delayed (the 
Carrington maps are created over ~27 days). Emerging regions and considerable 
changes in the global magnetic field may lead to errors in the resulting solar 
wind properties. However, these models provide a remarkable insight into the 
heliospheric morphology of the solar wind. The WSA model is currently used 
at the NOAA Space Weather Prediction Center and the forecast are available 
onlin^ The author suggests that these models could be aided by EUV or X-ray 
based CH maps, as it would increase the accuracy of the size and location of open 
magnetic fields, and thus may increase the accuracy of the model-based forecasts. 



A result of a recent solar wind forecast is shown in Figure 6.4 



^ http://www.swpc.noaa.gov/ws/ 
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Figure 6.4: The WSA solar wind velocity forecast. ACE solar wind speed data 
is shown with solid black line, while the colored dots show the WSA model pre- 
dictions. The red, green, and blue symbols correspond to predictions based on 
magnetograms obtained from Wilcox (WSO), Mount Wilson (MWO) and SOLIS 
(NSO) solar observatories, respectively. The uncertainty estimates are shown with 
vertical bars. 



A very exciting new development is the HELION project, which will provide 



access to data produced by different domains involved in solar and heliospheric 



http: / /www.helio-vo.eu/ aboutus/helio -consortium. php 
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6.2 Future work and outstanding questions 



physics. It will also allow the use of solar and heliospheric models in order to study 
solar events and phenomena in solar observations and link them to heliospheric 
effects, such as geomagnetic disturbances occurring on planets. The algorithms 
presented in this thesis will all be available through the HELIO development 
in the future. The CH Lambert maps and data products are already in the 
progress of being integrated into the HELIO database, while the HSSW profiles 
corresponding to CHs will also be publicly available in the coming months. These 



products will also be available on |www.SolarMonitor.org 
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